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ABSTRACT :  The  effect  of  forests  on  low  frequency  military  noise  propagation  is  unknown.  As  part  of  a  joint 
project,  ERDC-CERL  and  ERDC  CRREL  conducted  measurements  at  the  Lone  Star  Army  Ammunition  Plant 
located  in  Texarkana,  Texas,  to  investigate  these  effects.  In  this  report,  the  short-range  measurements  conducted  by 
ERDC-CRREL  are  documented.  Blast  noise  waveforms  produced  by  C4  explosions  at  distances  from  30  to  567  m 
were  recorded  and  are  presented  in  this  report.  In  all,  42  different  explosions  were  recorded,  producing  314  high 
quality  pressure  waveforms  for  analysis.  Additional  reports  documenting  the  long-range  measurements  and 
analyzing  the  recorded  data  are  in  preparation. 
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PREFACE 

This  report  was  prepared  by  Stephen  N.  Decato,  Physical  Science  Techni¬ 
cian,  Donald  G.  Albert,  Research  Geophysicist,  Frank  E.  Perron  Jr.,  Physical  Sci¬ 
ence  Technician,  and  David  L.  Carbee,  Physical  Science  Technician,  Geophysi¬ 
cal  Sciences  Branch,  U.S.  Army  Engineer  Research  and  Development  Center, 
Cold  Regions  Research  and  Engineering  Laboratory. 

In  recent  years  the  increase  of  civilian  populations  living  on  the  outskirts  of 
U.S.  Army  installations  has  led  to  more  frequent  noise  complaints.  While  the 
Army  normally  uses  forest  management  techniques  to  ensure  the  health  of  the 
forests  on  its  installations,  a  question  that  has  arisen  is  whether  these  forest  prac¬ 
tices,  including  tree  harvesting,  might  adversely  impact  the  noise  problem  by  re¬ 
ducing  the  noise  absorption  characteristics  of  the  forested  land.  Underlying  this 
question  is  the  assumption  that  forests  absorb  noise.  While  this  is  certainly  true  at 
high  frequencies  (in  the  kHz  band),  little  work  has  been  done  at  the  lower  fre¬ 
quencies  (below  100  Hz)  often  produced  by  Army  demolition  and  training  ac¬ 
tivities.  As  a  result  of  this  lack  of  knowledge,  ERDC,  under  the  direction  of  Dr. 
Larry  Pater,  ERDC-CERL,  is  embarking  on  a  study  of  the  low  frequency  noise 
attenuation  by  forests.  As  part  of  this  study,  a  joint  measurement  program  in¬ 
volving  two  ERDC  Laboratories,  CERL  and  CRREL,  was  conducted  at  the  Lone 
Star  Army  Ammunition  Plant  in  Texarkana,  Texas.  The  purpose  of  this  report  is 
to  document  the  CRREL  measurements.  Additional  reports  will  document  the 
CERL  measurements  and  analyze  the  experimental  data  to  determine  the  effect  of 
the  forest  on  low  frequency  military  blast  noise. 

Field  personnel  during  this  test  included  Don  Albert,  Steve  Decato,  and 
Frank  Perron  from  CRREL,  and  Larry  Pater,  Mike  White,  Richard  Racioppi,  Jeff 
Mifflin,  George  Swenson,  Brent  Miller,  and  Achal  Modi  from  CERL.  The  au¬ 
thors  especially  thank  Dave  Self,  Jesse  Stewart,  and  all  the  other  Lone  Star  Army 
Ammunition  Plant  personnel  for  their  assistance  during  these  measurements. 

They  thank  Tom  Vorac,  U.S.  AMC  Installations  and  Security  Activity  (now  with 
the  U.S.  Army  Environmental  Center,  Natural  Resources),  and  Dr.  Larry  Pater, 
U.S.A  ERDC-CERL,  for  funding  this  research.  Finally,  they  thank  Dr.  Pater, 
Project  Leader  for  Military  Noise  Management,  for  the  invitation  to  participate  in 
this  project. 

This  report  was  prepared  under  the  general  supervision  of  Dr.  Richard 
Detsch,  Chief,  Geophysical  Sciences  Branch;  and  James  Wuebben,  Acting 
Director,  CRREL. 

The  Commander  of  the  Engineer  Research  and  Development  Center  is  COL 
James  R.  Rowan,  EN.  The  Director  is  Dr.  James  R.  Houston. 
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1  INTRODUCTION 

This  report  documents  the  measurements  and  presents  the  signals  recorded 
during  the  Blast  Noise  Propagation  Through  Forest  Test  conducted  by  the  U.S. 
Army  ERDC  Construction  Engineering  Research  Laboratory  and  the  U.S.  Army 
ERDC  Cold  Regions  Research  and  Engineering  Laboratory.  The  measurements 
were  conducted  at  the  Lone  Star  Army  Ammunition  Plant  in  Texarkana,  Texas, 
from  23  July  through  25  July  2002  to  determine  the  properties  of  blast  noise 
wave  propagation  through  a  forested  area.  This  report  contains  the  recordings 
made  by  CRREL;  the  CERL  recordings  will  be  published  in  a  separate  report. 
Use  of  this  and  other  data  sets  will  assist  the  Army  in  determining  methods  of 
improving  blast  noise  attenuation  in  populated  areas  and  enhance  forest  man¬ 
agement  techniques. 

The  blast  waves  were  produced  by  detonating  charges  of  Ml  12  C4  explosive 
at  various  heights  and  locations.  Sledge  hammer  blows  on  an  aluminum  plate  and 
.45  caliber  blank  rounds  were  also  recorded  to  provide  additional  ground  and  air 
characterization.  Pressure  sensors,  microphones,  and  geophones  were  configured 
in  a  straight-line  array  on  the  ground  and  1.5  m  above  the  ground  surface. 

The  intent  of  this  report  is  to  provide  documentation  of  the  CRREL 
measurements  that  will  accompany  the  transmittal  of  the  recorded  data  on  CD- 
ROM  to  the  test  participants  and  sponsor.  This  report  includes  a  list  of  the 
records  obtained  during  the  tests,  documentation  of  the  sensor  array,  a  table  of 
sensors  used,  environmental  characterization,  and  a  description  of  the  data 
reduction.  Plots  of  all  the  signatures  are  provided,  and  the  signatures  themselves 
are  included  on  a  CD-ROM  in  ASCII  format.  An  analysis  report  of  the  signatures 
shown  here  will  be  published  separately. 
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2  OVERVIEW  OF  TEST 

The  Lone  Star  Army  Ammunition  plant  routinely  detonates  obsolete  muni¬ 
tions  as  a  method  of  disposal.  On  a  typical  day,  as  many  as  48  shots,  each  con¬ 
sisting  of  about  95  lb  (43  kg)  of  explosives,  may  be  detonated.  Although  these 
charges  are  buried  10  to  15  ft  (3  to  5  m)  deep  in  a  sandy/silty  soil,  the  shots  al¬ 
ways  “blow  out,”  sending  debris  high  into  the  air  producing  a  loud  blast  wave.  In 
recent  years,  private  homes  are  being  built  close  to  the  AAP  and  other  Army  fa¬ 
cilities’  boundary  lines.  This  has  resulted  in  concerns  about  noise  complaints 
from  new  residents.  As  the  current  demolition  area  at  the  AAP  is  surrounded  by 
forest,  a  scientifically  unanswered  question  is  what  is  the  effect  of  this  forest  on 
the  noise  produced  by  the  explosions?  If  the  trees  are  cut  or  thinned  as  a  part  of 
normal  forest  management  plan,  will  the  noise  complaints  increase?  To  answer 
this  question,  the  measurements  in  this  report  were  conducted. 

The  measurements  were  conducted  at  the  (new)  Demolition  Range  at  the 
Lone  Star  AAP.  Figure  1  shows  an  aerial  photograph  of  the  test  site  with  the  shot 
locations  (marked  TC  =  Test  Charge)  and  sensor  locations.  The  CRREL  sensor 
array  was  located  in  the  forest  south  of  TC  3.  Figure  2  shows  a  schematic  of  the 
CRREL  sensor  locations  and  a  typical  test  charge  (C-4)  shot.  Details  of  the  sen¬ 
sor  array  are  given  in  Section  3  of  this  report. 

Data  Recording  and  Reduction 

A  Geometries  NZ  digital  seismograph  (Fig.  3)  was  used  to  record  the  signa¬ 
ture  data  from  the  sensor  arrays.  We  used  an  8-kHz  sampling  rate  (0.125  ms)  per 
channel,  providing  a  3-dB  recording  bandwidth  of  1.75  Hz  to  3.3  kHz.  Twenty- 
four  sensor  channels  were  recorded  for  each  test  charge. 

Because  of  a  malfunction  with  the  CRREL  blaster  box,  the  actual  shot  instant 
could  not  be  recorded  and  the  recorder  could  not  be  automatically  started  as  we 
normally  do.  Instead,  the  recordings  were  started  manually  and  the  shot  instant 
was  determined  afterwards  from  the  blast  wave  arrival  times.  This  correction 
procedure  is  discussed  in  Section  6. 

The  seismograph  produces  binary  data  files  in  an  industry  standard  format 
called  SEG-2.  These  binary  files  contain  the  raw  voltage  output  of  the  sensors, 
but  contain  no  other  information.  We  converted  these  binary  files  into  ASCII 
files  using  the  MATLAB  program  listed  in  Appendix  B.  The  ASCII  files  have 
the  correct  physical  units,  time  information,  and  sensor  location  details.  The  sig¬ 
nature  plots  presented  in  Appendix  A  were  produced  by  the  MATLAB  program 
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listed  in  Appendix  C  from  the  ASCII  data  files.  Additional  data  analysis  will  be 
performed  using  only  the  ASCII  data  files. 


-600m 


-400m 


-200m 


200m 


Figure  1.  View  of  the  Lone  Star  AAP  test  range,  with  approximate  locations 
of  test  charge  and  sensor  locations. 
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Figure  2.  Test  Charge  Site  2  (TC  2)  and  a  sketch  of  the  test  array. 
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Figure  3.  Geometries  NZ  Digital  Seismograph,  set  up  inside  a 
vehicle  parked  behind  a  protection  soil  berm. 

List  of  Signature  Records 

The  measurements  were  conducted  without  interfering  with  the  normal 
demolition  operations.  This  meant  that  the  forest  testing  was  delayed  and  person¬ 
nel  were  evacuated  from  the  test  site  when  demolition  shots  were  fired.  While  we 
did  take  advantage  of  this  opportunity  to  record  some  of  these  demolition  shots 
using  a  separate  small  sensor  array  located  in  the  safety  zone,  those  recordings 
are  not  of  client  interest  to  this  study  and  are  not  included  in  this  report.  They  do, 
however,  affect  the  numbering  of  the  records  and  files  recorded  by  the  digital 
seismograph. 

The  following  list  contains  the  types  of  data  recorded  for  the  Blast  Noise 
Propagation  Test  and  is  included  on  the  CD.  NZ  file  names  like  “l.dat”  through 
“199.dat”  have  been  converted  to  multichannel  ASCII  files  with  the  correct 
physical  units.  File  names  like  “TX02002.asc”  through  “TX02046.asc”  are  files 
numbered  so  that  they  correspond  to  CERL’s  C4  shot  numbers.  Munition  demo¬ 
lition  records  have  been  omitted  from  this  report.  Table  1  contains  the  list  of  sig¬ 
nature  records.  Figures  4  and  5  show  our  C4  charges  and  .45-calibre  pistol. 
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Table  1.  Signature  records. 


1  dat-  19  dat 

.45  cal.  bLank 

20.dat  -  27.dat 

Microphone  calibration 

28.dat-  3S.dat 

Demolition  of  munitions 

39.dat-44.dat 

C4 

45.dat 

Noise  Record 

46.dat-  56.dat 

Demolition  of  munitions 

57.dat-58.dat 

Noise  Record 

59.dat 

Demolition  of  munitions 

60.dat  -68.dat 

C4 

69.dat 

Noise  Record 

70.dat-  74.dat 

C4 

75.dat-81.dat 

Demolition  of  munitions 

82.dat-91.dat 

Calibration 

92.dat 

Noise  Record 

93.dat-113.dat 

C4 

1 14.dat  -  199  dat 

Calibration 

Specifications  for  Sources  used  in  Measurements 
C4 


.45  cal  blank 


V2  brick 

0.63  lb  (0.28  kg) 

1  brick 

1.25  lb  (0.57  kg) 

2  bricks 

2.50  lb  (1.23  kg) 

4  bricks 

5.00  lb  (2.27  kg) 

5  bricks 

6.25  lb  (2.84  kg) 

8  bricks 

10.00  lb  (4.55  kg) 

10  bricks 

12.50  lb  (5.67  kg) 

Overall  Length 

0.89  in.  (22.6  mm) 

Case  Length 

0.89  in.  (22.6  mm) 

Cartridge  Weight  7.6  g 

Case  Weight 

5.8  g 

Powder  Weight 

31  g 

Powder  Type 

Black,  FFFFG,  CCI 

Primer  Type 

Magnum 

Manufacture 

Custom  made 

Overall  Length 

34  in.  (86.7  cm) 

Overall  Weight 

11.68  lb  (5.30  kg) 

Head  Weight 

10.0  lb  (4.55  kg) 

Strike  Plate  Material 
Strike  Plate  Area 
Strike  Plate  Thickness 
Strike  Plate  Weight 


606 1T6  Aluminum 
103.1  in.2  (668  cm2) 
1  in.  (2.54  cm) 
10.36  lb.  (4.70  kg) 
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Figure  4.  .45-caliber  calibration  pistol  firing  blanks  at  a  1-m  height. 


a.  C4  being  formed  into  a  5-lb  (2.25-kg)  spherical  charge. 
Figure  5.  Making  the  C4  charges. 
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b.  A  1.25-lb  (0.68-kg)  brick  of  C4. 


c.  C4  spherical  charge  being  suspended  at  TC  3. 
Figure  5  (con't).  Making  the  C4  charges. 
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3  CRREL  SENSOR  CONFIGURATION 

The  sensor  configuration  installed  by  CRREL  consisted  of  a  linear  array  with 
multi-sensor  stations.  The  array  was  located  in  the  forest  between  TC  3  and  TC  4. 
It  extended  from  30  to  151.3  m  magnetic  south  of  TC  3.  The  goal  of  this  sensor 
array  was  to  obtain  a  detailed  view  of  the  blast  waves  as  they  propagated  into  the 
forest. 

CERL  also  installed  and  recorded  signatures  from  their  own  independent 
sensor  array.  This  array  was  located  both  inside  and  outside  of  the  forest  on  a 
true  north-south  line.  This  array  was  designed  to  obtain  a  complete  overview  of 
the  blast  wave  propagation  at  the  site  as  well  as  details  of  the  blast  wave 
generation  by  the  source. 

Table  2  lists  the  distances  between  the  CRREL  and  CERL  locations.  This 
difference  was  due  to  a  3.5  degree  east  deviation  between  true  north  and 
magnetic  north. 

Figures  6  through  1 1  show  the  locations  of  the  of  the  CRREL  instrument 
stations.  Table  3  details  the  sensors  used  in  the  CRREL  Blast  Wave  Propagation 
Array.  Table  4  notes  any  sensor  changes  made  to  that  array,  during  the 
measurements. 

Figure  12  shows  the  CRREL  array  used  to  record  the  large  munitions 
demolition  shots.  Further  discussion  of  this  array  is  omitted  from  this  report. 


Table  2.  Distances  between  the  CRREL  and  CERL  locations. 


Location 

Distance  (m) 

30  meters 

1.83 

60  meters 

3.66 

90  meters 

5.49 

120  meters 

7.32 

150  meters 

9.15 
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Figure  6.  30  meter  sensor  locations. 


Seismic  and  Acoustic  Signature  Measurements  Through  Forest 


11 


TCI 

(-414.7  m) 


North 


North  (M) 


f 


Array  #  1 


m/ 


TC2 

(-60.5  tn) 


,/t 


TC3 
(0.0  m) 


—  30  m 
-""60  m 

—  90  m 
""120  m 

150  m  &  151.3  m 


TC4 
(311  m) 


•  =  CRIREL  Sensor  Sites 
O  =  CERL  Sites 


Figure  7.  60  meter  sensor  locations. 
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Figure  8.  90  meter  sensor  locations, 
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Figure  9.  120  meter  sensor  locations. 
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Figure  10. 150  meter  sensor  locations. 
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Figure  11.  153.1  meter  sensor  locations. 
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Figure  12.  CRREL  munitions  demolition  array.  This  array 
was  located  approximately  600  to  900  m  from  the  demo¬ 
lition  charges,  with  Station  4  closest  to  the  charges. 
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Table  3.  Wave  propagation  sensor  array. 


Array  Cable  1 


Location 

Channel 

Number 

Sensor 

Serial 

Number 

Height 

(m) 

Calibrated  Sensitivity 

30  Meter 

1 

Radial  geophone 

H13 

0.0 

32.2  V/m/s 

30  Meter 

2 

Vertical  geophone 

V23 

0.0 

32.2  V/m/s 

30  Meter 

3 

PCB  102A07 

15972 

0.0 

1.28E-02  mV/Pa 

30  Meter 

4 

PCB  102A07 

15973 

1.5 

1.24E-02  mV/Pa 

60  Meter 

5 

Radial  geophone 

H2 

0.0 

32.2  V/m/s 

60  Meter 

6 

Vertical  geophone 

VII 

0.0 

32.2  V/m/s 

60  Meter 

7 

PCB  102A07 

15971 

0.0 

1.30E-02  mV/Pa 

60  Meter 

8 

PCB  106B50 

6693 

1.5 

6.74E-02  mV/Pa 

90  Meter 

9 

Radial  geophone 

H18 

0.0 

32.2  V/m/s 

90  Meter 

10 

Vertical  geophone 

V4 

0.0 

32.2  V/m/s 

90  Meter 

11 

PCB  106B50 

3259 

0.0 

8.17E-02  mV/Pa 

90  Meter 

12 

PCB  106B50 

6695 

1.5 

7.16E-02  mV/Pa  (uncal) 

Array  Cable  2 


Location 

Channel 

Number 

Sensor  Number 

Sensor  S/N 

Height 

(m) 

Calibrated  Sensitivity 

120  Meter 

13 

Radial  geophone 

H28 

0.0 

32.2  V/m/s 

120  Meter 

14 

Transv.  geophone 

H19 

0.0 

32.2  V/m/s 

120  Meter 

15 

Vertical  geophone 

V6 

0.0 

32.2  V/m/s 

120  Meter 

16 

B&K  4938  -  3 

2239255 

0.0 

1.15  E-3  mV/Pa 

150  Meter 

17 

Radial  geophone 

H20 

0.0 

32.2  V/m/s 

150  Meter 

18 

Transv.  geophone 

H21 

0.0 

32.2  V/m/s 

150  Meter 

19 

Vertical  geophone 

V28 

0.0 

32.2  V/m/s 

150  Meter 

20 

B&K  4938  -  2 

2239254 

0.0 

1.32  E-3  mV/Pa 

151.3  Meter 

21 

Radial  geophone 

H10 

1.5 

32.2  V/m/s 

151.3  Meter 

22 

Transv.  geophone 

H9 

1.5 

32.2  V/m/s 

151.3  Meter 

23 

Vertical  geophone 

V? 

1.5 

32.2  V/m/s 

151.3  Meter 

24 

B&K  4938  -  1 

2239253 

1.5 

1 .20  E-3  mV/Pa 

All  geophones  are  Mark  Products  L-15,  4.5  Hz 
All  B&K  microphones  had  1:1000  attenuation. 


18 


ERDC/CRREL  TR-05-10 


Table  4.  Changes  to  wave  propagation  sensor  array. 


The  PCB  sensor  at  the  60  m  location  failed  and  was  replaced  after  CRREL  binary 
file  number  44,  or  CERL  shot  number  6. 


Location 

Channel 

Sensor  Number 

Sensor  S/N 

Height  |  Calibrated  Sensitivity 

60  Meter 

7 

PCB  106B50 

3260 

0.0  7.16E-02  mV/Pa 

60  Meter 

7 

PCB  102A07 

15971 

0.0  1.30E-02  mV/Pa 

A  =  Binary  records  45  -  199 
B  =  Binary  records  1  -  44 


The  B&K  4938  microphones  at  locations  120  m  and  150  m  were  replaced  with 
with  B&K  4165  microphones  for  blank  pistol  shots  5-19,  ascii 
files  105  - 119,  plus  the  associated  calibration  files. 


Location 

Channel 

Number 

Sensor  Number 

Sensor  S/N 

Height 

(m) 

Calibrated  Sensitivity 

120  Meter 

16  A 

B&K  4938  -  3 

2239255 

0.0 

1.15  E-3  mV/Pa 

120  Meter 

16  B 

B&K  4165 -5 

1881043 

0.0 

0.75  E-3  mV/Pa 

150  Meter 

20  A 

B&K  4938  -  2 

2239254 

0.0 

1.32  E-3  mV/Pa 

1 50  Meter 

20  B 

B&K  4165 -6 

1857571 

0.0 

0.73  E-3  mV/Pa 

A  =  Binary  records  39-199 
B  =  Binary  records  5  -  27 

The  4165  microphones  had  too  low  of  a  calibration  value,  so  may  not  be  reliable. 
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4  SENSOR  CALIBRATION 

Tables  3  and  4  list  the  sensitivities  used  to  convert  the  voltages  recorded  by 
the  seismograph  to  physical  units.  These  two  tables  also  contain  the  calibration 
values  for  all  sensors  used.  The  sensitivities  were  determined  as  follows: 

■  Geophones :  The  sensitivity  provided  by  the  manufacturer  was  used.  From 
past  experience  at  CRREL,  these  sensors  are  usually  accurate  to  about  10%, 
but  they  can  occasionally  vary  more  than  this  range  because  of  differences  in 
coupling  to  the  ground. 

•  Pressure  Sensors  (solid  state):  The  PCB  pressure  sensors  were  calibrated  in 
situ  using  a  Trig-Tek  Model  402H  calibrator  (Fig.  13).  This  model  calibrator 
provides  sound  pressure  levels  of  120  to  160  dB  in  10-dB  increments.  We 
used  levels  of  140,  150,  and  160  dB  for  these  calibrations  to  match  the 
expected  blast  noise  levels.  It  has  selectable  internal  low  frequency  outputs 
of  62.5,  125,  250,  500  and  1000  Hz.  We  relied  mostly  on  the  lower 
frequencies,  62.5-250  Hz,  for  these  calibrations. 

•  Microphones  (  Va  and  Vi  inch  B&K):  The  B&K  %  inch  Model  4938 
microphones  were  also  calibrated  using  the  Trig-Tek  Model  402H  calibrator 
at  140  and  150  dB. 


Figure  13.  Typical  sensor  calibration  using  the  Trig-Tek 
Model  402 H  calibrator. 
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The  B&K  !4  inch  Model  4165  microphones  were  calibrated  using  a  GenRad 
Type  1562-A  sound  level  calibrator.  This  type  of  calibrator  provides  a  sound 
pressure  level  of  1 14  dB.  It  has  selectable  internal  frequency  outputs  of  125,  250, 
500,  1000  and  2000  Hz. 
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5  FOREST  ENVIRONMENT 

Figures  13  through  19  show  typical  views  of  the  conditions  within  the  forest. 
The  understory  vegetation  limited  visibility  to  about  30-60  m  within  the  forest. 
Details  of  the  tree  species  and  density  are  to  be  provided  by  forest  measurements 
already  conducted  as  part  of  the  land  management  at  the  Lone  Star  AAP. 

Ground  conditions  usually  consisted  of  a  layer  of  needles  a  few  inches  thick, 
as  seen  in  the  photographs.  The  weather  during  the  tests  was  very  hot  and  humid, 
with  air  temperatures  typically  near  100°F  (38°C). 


Figure  14.  Vegetation  at  30  meter  site. 
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Figure  15.  Vegetation  at  60  m  site. 


Figure  16.  Vegetation  at  90  m  site 
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Figure  17.  Vegetation  at  120  m  site. 


Figure  18.  Vegetation  at  150  m  site. 
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Figure  19.  Vegetation  at  151.3  m  site 
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6  TEXAS  2002  TIME  BREAK  CORRECTIONS  AND 
ACOUSTIC  VELOCITIES 

Because  of  a  malfunction  with  the  blaster  box,  good  time  breaks  were  not 
obtained.  The  seismograph  recordings  were  started  manually  during  the 
countdown,  so  the  blast  wave  arrivals  occur  at  different  times  within  the  binary 
file. 


Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 


Before  writing  the  ASCII  files,  the  true  time  break  was  found  by  finding  the 
blast  wave  arrival  times  and  the  correct  zero  time  from  the  acoustic  velocity. 
These  corrections  are  listed  below.  In  the  ASCII  files,  the  blast  wave  arrivals 
occur  at  the  true  arrival  time. 

In  this  listing,  the  shot  number  is  the  original  CRREL  binary  file  (xx.dat) 
recorded  by  the  NZ  seismograph.  For  C4  shots,  the  ASCII  file  name  with  the 
time  break  corrections  applied  is  tx020xx.asc,  where  xx  is  the  CERL  number 
given  below.  For  pistol  shots,  the 

•  ASCII  file  name  is  tx021xx.asc,  where  xx  is  the  CRREL  shot  number. 

•  (For  shot  40,  the  binary  file  is  40.dat  and  the  ascii  file  is  tx2002.asc.) 

•  (For  shot  1,  the  binary  file  is  l.dat  and  the  ascii  file  is  tx2101.asc.) 

The  acoustic  velocities  listed  below  for  the  pistol  shots  are  higher  than  the 
true  velocities  because  of  the  offset  of  the  pistol  from  the  line  of  the  sensors.  As 
the  pistol  was  only  2  m  away  from  the  nearest  sensor,  this  did  not  affect  the  time 
break  correction. 


The  records  are  sorted  by  shot  location  below. 


60  TC2  CERL  7  Velocity  =  349.7 
62  TC2  CERL  9  Velocity  =  348.0 
66  TC2  CERL  13  Velocity  =  348.2 

70  TC2  CERL  3B  Velocity  =  347.4 

71  TC2  CERL  4B  Velocity  =  347.3 

72  TC2  CERL  5B  Velocity  =  348.2 

73  TC2  CERL  6B  Velocity  =  348.7 

74  TC2  CERL  7B  Velocity  =  348.6 

94  TC2  CERL  17  Velocity  =  349.6 

97  TC2  CERL  21  Velocity  =  350.7 

101  TC2  CERL  25  Velocity  =  350. 

105  TC2  CERL  29  Velocity  =  351. 

108  TC2  CERL  32  Velocity  =  352. 


m/s  TB  corr  =  -0.0465  7-24-02 
m/s  TB  corr  =  -2.5717  7-24-02 
m/s  TB  corr  =  0.0052  7-24-02 
m/s  TB  corr  =  -1.2592  7-25-02 
m/s  TB  corr  =  -1.0796  7-25-02 
m/s  TB  corr  =  -1.3115  7-25-02 
m/s  TB  corr  =  -1.2592  7-25-02 
m/s  TB  corr  =  -1.0384  7-25-02 
m/s  TB  corr  =  -0.3061  7-25-02 
m/s  TB  corr  =  -0.2621  7-25-02 
5  m/s  TB  corr  =  -0.5556  7-25-02 

0  m/s  TB  corr  =  -0.6646  7-25-02 

3  m/s  TB  corr  =  -0.8483  7-25-02 
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Shot  109  TC2  CERL  33  Velocity  =  352.7  m/s  TB  corr  =  -1.2742  7-25-02 

Shot  110  TC2  CERL  34  Velocity  =  351.8  m/s  TB  corr  =  -1.5879  7-25-02 

Shot  111  TC2  CERL  35  Velocity  =  352.4  m/s  TB  corr  =  -1.0272  7-25-02 

Shot  112  TC2  CERL  36  Velocity  =  351.9  m/s  TB  corr  =  -0.9628  7-25-02 

Shot  113  TC2  CERL  37  Velocity  =  351.5  m/s  TB  corr  =  -1.2390  7-25-02 

Shot  61  TCI  CERL  8  Velocity  =  350.2  m/s  TB  corr  =  1.2226  7-24-02 
Shot  65  TCI  CERL  12  Velocity  =  348.6  m/s  TB  corr  =  0.2675  7-24-02 
Shot  93  TCI  CERL  16  Velocity  =  352.9  m/s  TB  corr  =  -0.4247  7-25-02 

Shot  96  TCI  CERL  20  Velocity  =  348.7  m/s  TB  corr  =  -0.2756  7-25-02 

Shot  100  TCI  CERL  24  Velocity  =  351.5  m/s  TB  corr  =  -0.6963  7-25-02 

Shot  104  TCI  CERL  28  Velocity  =  351.0  m/s  TB  corr  =  -0.8121  7-25-02 

Shot  40  TC3  CERL  2  Velocity  =  350.4  m/s  TB  corr  =  -0.2134  7-24-02 

Shot  41  TC3  CERL  3  Velocity  =  349.0  m/s  TB  corr  =  -0.2179  7-24-02 

Shot  42  TC3  CERL  4  Velocity  =  350.5  m/s  TB  corr  =  -0.2174  7-24-02 

Shot  67  TC3  CERL  14  Velocity  =  348.7  m/s  TB  corr  =  0.0880  7-24-02 

Shot  98  TC3  CERL  22  Velocity  =  350.6  m/s  TB  corr  =  -1.4086  7-25-02 
Shot  102  TC3  CERL  26  Velocity  =  350.4  m/s  TB  corr  =  -0.7076  7-25-02 

Shot  106  TC3  CERL  30  Velocity  =  350.9  m/s  TB  corr  =  -0.4905  7-25-02 

Shot  64  TC4  CERL  11  Velocity  =  353.1  m/s  TB  corr  =  -0.8379  7-24-02 

Shot  68  TC4  CERL  15  Velocity  =  353.3  m/s  TB  corr  =  -1.0028  7-24-02 

Shot  95  TC4  CERL  19  Velocity  =  355.0  m/s  TB  corr  =  -1.1580  7-25-02 

Shot  99  TC4  CERL  23  Velocity  =  355.2  m/s  TB  corr  =  -1.2413  7-25-02 

Shot  103  TC4  CERL  27  Velocity  =  355.4  m/s  TB  corr  =  -1.5833  7-25-02 

Shot  107  TC4  CERL  31  Velocity  =  355.7  m/s  TB  corr  =  -2.0720  7-25-02 

Blank  Pistol  Shots 

(The  velocity  is  an  apparent  phase  velocity,  not  the  true  acoustic  velocity,  because  the  source 
was  offset  from  the  sensor  line.) 

Shot  1  30m  Velocity  =  362.3  m/s  TB  corr  =  0.8091  7-23-02 

Shot  2  30m  Velocity  =  362.5  m/s  TB  corr  =  0.5928  7-23-02 

Shot  3  30m  Velocity  =  362.8  m/s  TB  corr  =  0.9599  7-23-02 

Shot  4  30m  Velocity  =  363.9  m/s  TB  corr  =  1.2268  7-23-02 

Shot  15  30m  Velocity  =  362.3  m/s  TB  corr  =  0.9330  7-24-02 

Shot  16  30m  Velocity  =  362.0  m/s  TB  corr  =  1.2151  7-24-02 

Shot  17  30m  Velocity  =  362.0  m/s  TB  corr  =  1.3954  7-24-02 

Shot  18  30m  Velocity  =  362.0  m/s  TB  corr  =  1.1648  7-24-02 

Shot  19  30m  Velocity  =  362.3  m/s  TB  corr  =  1.1863  7-24-02 
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Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 

Shot 


5  90m  Velocity  = 

6  90m  Velocity  = 

7  90m  Velocity  = 

8  90m  Velocity  = 

9  90m  Velocity  = 

10  60m  Velocity  = 

11  60m  Velocity  = 

12  60m  Velocity  = 

13  60m  Velocity  = 

14  60m  Velocity  = 


356.9  m/s  TB  corr  = 
359.0  m/s  TB  corr  = 
358.2  m/s  TB  corr  = 
358.5  m/s  TB  corr  = 
359.0  m/s  TB  corr  = 

364.2  m/s  TB  corr  = 
364.5  m/s  TB  corr  = 
364.5  m/s  TB  corr  = 
364.5  m/s  TB  corr  = 
364.5  m/s  TB  corr  = 


1.4087  7-24-02 
1.4700  7-24-02 
1.0701  7-24-02 
1.0898  7-24-02 
1.2376  7-24-02 

1.2334  7-24-02 
1.2090  7-24-02 
1.2828  7-24-02 
1.0293  7-24-02 
1.3029  7-24-02 
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7  SIGNATURE  RECORDINGS  AND  DATA  QUALITY 

This  section  provides  plots  of  all  of  the  C4  test  charges  and  all  of  the  .45 
caliber  pistol  shots  recorded  by  the  CRREL  sensor  array.  The  plots  were  pro¬ 
duced  using  the  ASCII  data  files  and  the  MATLAB  plot  program  listed  in 
Appendix  C. 

Table  5  lists  the  events  recorded  during  the  tests.  As  explained  in  Appendix 
B,  the  binary  data  file  names  were  changed  to  agree  with  the  CERL  Test  Charge 
numbering  system.  Also,  to  avoid  duplicate  file  names,  the  .45  caliber  shot  files 
were  renumbered  in  the  range  101-119. 

Table  6  gives  an  assessment  of  the  data  quality  of  each  trace,  based  on  the 
plots  in  this  section.  The  table  lists  “good”  data  suitable  for  additional  analysis. 
“Bad”  data  should  not  be  used,  at  least  not  without  further  processing.  The  data 
might  be  bad  because  of  sensor  failure  or  too  high  a  noise  level. 

Nine  pressure  sensor  channels  were  recorded  for  each  of  the  C4  shots.  For 
the  42  C4  shots  detonated  during  the  test,  five  shots  were  not  recorded  because  of 
a  failure  with  the  time  break  signal.  Four  other  channels  were  unusable,  while  15 
were  usable  but  noisy. 

Appendix  A  shows  each  individual  record  that  was  made. 


Seismic  and  Acoustic  Signature  Measurements  Through  Forest 


29 


Table  5.  Chronology  of  recorded  events. 

CRREL  CERL  ASCII 

File  Shot  File  Local  Source  Ht. 


Name 

# 

Name 

Date 

Time 

Source 

m  (ft/in.) 

Comments 

1.dat 

TX02101.asc 

07/23/02 

1545 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

2.dat 

TX02102.asc 

07/23/02 

1553 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

[  3.dat 

TX02103.asc 

07/23/02 

1555 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

4.dat 

TX02104.asc 

07/23/02 

1556 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

5.dat 

TX02105.asc 

07/24/02 

0749 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  west  of  90  m  site 

6.dat 

TX02106.asc 

07/24/02 

0753 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  west  of  90  m  site 

7.dat 

TX02107.asc 

07/24/02 

0754 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  west  of  90  m  site 

8.dat 

TX02108.asc 

07/24/02 

0755 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  west  of  90  m  site 

9.dat 

TX02109.asc 

07/24/02 

0756 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  west  of  90  m  site 

10.dat 

TX021  lO.asc 

07/24/02 

0807 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  60  m  site 

1 1  .dat 

TX02111.asc 

07/24/02 

0809 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  60  m  site 

12.dat 

TX021 12.asc 

07/24/02 

0810 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  60  m  site 

13.dat 

TX021 13.asc 

07/24/02 

0811 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  60  m  site 

14.dat 

TX021 14.asc 

07/24/02 

0811 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  60  m  site 

15.dat 

TX021 15.asc 

07/24/02 

0817 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

16.dat 

TX021 16.asc 

07/24/02 

0818 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

17.dat 

TX021 17.asc 

07/24/02 

0818 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

18.dat 

TX021 18.asc 

07/24/02 

0819 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

19.dat 

TX021 19.asc 

07/24/02 

0819 

.45  cal.  blank 

1.0  (3/3) 

Shot  2m  east  of  30  m  site 

20.dat 

07/24/02 

Calibration 

20.dat  thru  27.dat  are  calibration  files 

27.dat 

07/24/02 

0909 

Calibration 

28.dat 

07/24/02 

1009 

Munition  Demo 

28.dat  thru  38.dat  are  munition  demo  files 

38.dat 

07/24/02 

1014 

Munition  Demo 

39.dat 

1 

07/24/02 

1218 

C4  (1  brick) 

1.9  (6/3) 

TC2  (See  Note  1) 

40.dat 

2 

TX02002.asc 

07/24/02 

1235 

C4  (4  bricks) 

1.9  (6/3) 

TC2 

41  .dat 

3 

TX02003.asc 

07/24/02 

1334 

C4  (1  brick) 

3.8  (12/6) 

TC2 

42.dat 

4 

TX02004.asc 

07/24/02 

1345 

C4  (1  brick) 

3.8  (12/6) 

TC2 

43.dat 

5 

07/24/02 

1356 

C4  (1  brick) 

1.2  (4/0) 

TC2  (See  Note  1) 

44.dat 

6 

07/24/02 

1406 

C4  (1  brick) 

0.6  (2/0) 

TC2  (See  Note  1) 

45.dat 

07/24/02 

1515 

Noise  Record 

46.dat 

07/24/02 

Munition  Demo 

46.dat  thru  56.dat  are  munition  demo  files 

56.dat 

07/24/02 

Munition  Demo 

57.dat 

07/24/02 

Noise  Record 

58.dat 

07/24/02 

Noise  Record 

59.dat 

07/24/02 

Munition  Demo 

Note  1 .  CRREL  shotbox  sent  early  time  break.  No  ascii  file. 


30 


ERDC/CRREL  TR-05-10 


Table  5  (cont'd).  Recorded  events. 


CRREL  CERL  ASCII 

File  Shot  File  Local  Source  Ht. 


Name 

# 

Name 

Date 

Time 

Source 

m  (ft/in.) 

Comments 

60.dat 

7 

TX02007.asc 

07/24/02 

1550 

C4  (1  brick) 

0.3  (1/0) 

TC2 

61  .dat 

8 

TX02008.asc 

07/24/02 

1613 

C4  (1  brick) 

1.9  (6/3) 

TCI 

62.dat 

9 

TX02009.asc 

07/24/02 

1614 

C4  (1  brick) 

1.9  (6/3) 

TC2 

63.dat 

10 

07/24/02 

1615 

C4  (1  brick) 

1.9  (6/3) 

TC3  (See  Note  1) 

64.dat 

11 

TX02011.asc 

07/24/02 

1617 

C4  (1  brick) 

1.9  (6/3) 

TC4 

65.dat 

12 

TX02012.asc 

07/24/02 

1631 

C4  (1  brick) 

1.9  (6/3) 

TCI 

66.dat 

13 

TX02013.asc 

07/24/02 

1633 

C4  (1  brick) 

1.9  (6/3) 

TC2 

67.dat 

14 

TX02014.asc 

07/24/02 

1634 

C4  (1  brick) 

1.9  (6/3) 

TC3 

68.dat 

15 

TX02015.asc 

07/24/02 

1635 

C4  (1  brick) 

1.9  (6/3) 

TC4 

69.dat 

07/25/02 

0820 

Noise  Record 

70.dat 

43 

TX02043.asc 

07/25/02 

0859 

C4  (1  brick) 

3.8  (12/6) 

TC2 

Repeat  of  CERL  shot  #  3 

71  .dat 

44 

TX02044.asc 

07/25/02 

0910 

C4  (1  brick) 

3.8  (12/6) 

TC2 

Repeat  of  CERL  shot  #  4 

72.dat 

45 

TX02045.asc 

07/25/02 

0933 

C4  (1  brick) 

3.0  (10/0) 

TC2 

Repeat  of  CERL  shot  #  5 

73.dat 

46 

TX02046.asc 

07/25/02 

0942 

C4  (1  brick) 

1.2  (4/0) 

TC2 

Repeat  of  CERL  shot  #  6 

74.dat 

47 

TX02047.asc 

07/25/02 

0951 

C4  (1  brick) 

0.3  (1/0) 

TC2 

Repeat  of  CERL  shot  #  7 

75.dat 

07/25/02 

1005 

Munition  Demo 

75.dat  thru  81  .dat  are  munition  demo  files 

81.dat 

07/25/02 

Munition  Demo 

82.dat 

07/25/02 

1013 

Calibration 

82.dat  thru  91  .dat  are  munition  demo  files 

91  .dat 

07/25/02 

1019 

Calibration 

92.dat 

07/25/02 

1103 

Noise  Record 

93.dat 

16 

TX02016.asc 

07/25/02 

1106 

C4  (1  brick) 

1.9  (6/3) 

TCI 

94.dat 

17 

TX02017.asc 

07/25/02 

1107 

C4  (1  brick) 

1.9  (6/3) 

TC2 

18 

07/25/02 

1107 

C4  (1  brick) 

1.9  (6/3) 

TC3  (See  note  2) 

95.dat 

19 

TX02019.asc 

07/25/02 

1108 

C4  (1  brick) 

1.9  (6/3) 

TC4 

96.dat 

20 

TX02020.asc 

07/25/02 

1122 

C4  (5  pounds) 

1.9  (6/3) 

TCI  (See  note  3) 

97.dat 

21 

TX02021  .asc 

07/25/02 

1123 

C4  (5  pounds) 

1.9  (6/3) 

TC2 

98.dat 

22 

TX02022.asc 

07/25/02 

1124 

C4  (5  pounds) 

1.9  (6/3) 

TC3 

99.9at 

23 

TX02023.asc 

07/25/02 

1125 

C4  (5  pounds) 

1.9  (6/3) 

TC4 

100.dat 

24 

TX02024.asc 

07/25/02 

1143 

C4  (5  pounds) 

1.9  (6/3) 

TCI 

101.dat 

25 

TX02025.asc 

07/25/02 

1144 

C4  (5  pounds) 

1.9  (6/3) 

TC2 

102.dat 

26 

TX02026.asc 

07/25/02 

1145 

C4  (5  pounds) 

1.9  (6/3) 

TC3 

103.dat 

27 

TX02027.asc 

07/25/02 

1146 

C4  (5  pounds) 

1.9  (6/3) 

TC4 

104.dat 

28 

TX02028.asc 

07/25/02 

1202 

C4  (5  pounds) 

1.9  (6/3) 

TCI 

105.dat 

29 

TX02029.asc 

07/25/02 

1203 

C4  (5  pounds) 

1.9  (6/3) 

TC2 

106.dat 

30 

TX02030.asc 

07/25/02 

1204 

C4  (5  pounds) 

1.9  (6/3) 

TC3 

Note  2.  Time  break  failed.  No  binary  or  ascii  file. 


Note  3.  All  C-4  charges  detonated  to  create  files  96.dat  through  1 13.dat  were  done  using  spherical  charges,  not  preformed  bricks. 
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Table  5  (cont'd). 


CRREL  CERL  ASCII 

File  Shot  File  Local  Source  Ht. 


Name 

# 

Name 

Date 

Time 

Source 

m  (ft/in.) 

Comments 

107.dat 

31 

TX02031.asc 

07/25/02 

1204 

C4  (5  pounds) 

1.9  (6/3) 

TC4 

108.dat 

32 

TX02032.asc 

07/25/02 

1316 

C4  (5  pounds) 

1.9  (6/3) 

TC2 

109.dat 

33 

TX02033.asc 

07/25/02 

1325 

C4  (5  pounds) 

3.8  (12/6) 

TC2 

110.dat 

34 

TX02034.asc 

07/25/02 

1336 

C4  (5  pounds) 

3.8  (12/6) 

TC2 

1 1 1  .dat 

35 

TX02035.asc 

07/25/02 

1347 

C4  (5  pounds) 

3.0  (10/0) 

TC2 

112.dat 

36 

TX02036.asc 

07/25/02 

1356 

C4  (5  pounds) 

1.2  (4/0) 

TC2 

113.dat 

37 

TX02037.asc 

07/25/02 

1406 

C4  (5  pounds) 

0.3  (1/0) 

TC2 

114.dat 

07/25/02 

1543 

Calibration 

1 14.dat  thru  199.dat  are  calibration  files 

199.dat 

07/25/02 

1648 

Calibration 

Table  6.  Data  quality. 


TX  2002  Forest  Test 
Sensors  -  Channel  numbers 


C4  Shots 


CERL 


Rec  Num  -> 

1  2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

CRREL 

Rec  Num-> 

39  40 

41 

42 

43 

44 

60 

61 

62 

63 

64 

65 

66 

67 

68 

93 

94 

X 

95 

96 

97 

98 

99 

4  P-30-1 
3  P-30-0 
8  P-60-1 
7  P-60-0 
12  P-90-1 
1 1  P-90-0 
16  P-120-0 
20  P-150-0 
24  P-152T-1 


Shot  box  TB  failed  for  CERL  shots  1,  5,  6,  10.  No  ascii  files. 
TB  failure  for  CERL  shot  18,  no  binary  or  ascii  file. 
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Table  6  (cont'd). 

C4  Shots 


3B  4B  5B  6B  7B 


CERL 

Rec  Num  -> 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

43 

44 

45 

46 

47 

CRREL 

Rec  Num-> 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

70 

71 

72 

73 

74 

4  P-30-1 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

3  P-30-0 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

8  P-60-1 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

7  P-60-0 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

12  P-90-1 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

1 1  P-90-0 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

16  P-120-0 

O 

O 

O 

O 

N 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

20  P-150-0 

O 

O 

O 

O 

N 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

24  P-152T-1 

O 

O 

O 

O 

N 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O  =  Good 
X=  Bad 

N  =  Noisy  data  -  but  usable 
R  =  Ringing  sensor 
C  =  Crosstalk 


BLANK  PISTOL  SHOTS 


CERL 

Rec  Num  -> 
CRREL 
Rec  Num-> 

4  P-30-1 
3  P-30-0 
8  P-60-1 
7  P-60-0 
12  P-90-1 
11  P-90-0 
16  P-120-0 
20  P-150-0 
24  P-152T-1 


101  102  103  104  105  106  107  108  109  110  111  112  113  114  115  116  117  118  119 


1 

O 

O 

O 

N 

O 

N 

X 

X 

X 


10  11 
X  X 


12  13  14  15  16  17  18 


19 

O 

X 

O 

X 

O 

N 

X 

X 

X 


For  records  105-119,  channels  16,  20  used  B&K4165  microphones. 


O  =  Good 
X=  Bad 

N  =  Noisy  data  -  but  usable 
R  =  Ringing  sensor 
C  =  Crosstalk 
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APPENDIX  A:  SIGNAL  DATA  PLOTS 


The  following  signal  data  plots  show  the  response  of  the  24  sensors  recorded 
during  each  test.  The  plots  are  normalized  so  that  each  channel  is  the  same  size. 
The  maximum  amplitude  is  shown  in  the  label  for  each  channel.  Figure  A1 
describes  the  plot  layout,  labels  and  descriptions  used  for  each  record  made 
during  the  test.  The  plots  are  normalized  so  that  each  channel  is  the  same  size. 
The  maximum  amplitude  is  shown  in  the  label  for  each  channel. 


i  Project  j 

File  |  |  Test  Charge  ! 

i  i 

!  Source  and 

i 

!  Distance  from  source  to  i 

]  Name  \ 

I  Name  [  Location 

|  Height  | 

all  sensor  locations 

i  +  *  , - - - ,  i 

TX02  —  TX02002.asc  —  TC3  2.3  kg  C-4  1.9  magi  —  30-152  m 


Figure  A1.  Description  of  the  signal  data  plots. 


TX02002.asc  —  TC3  2.3  kg  C-4  1.9  m  agl  —  30-152  m 


Time,  sec  Time,  sec  Time,  sec 


TX02003.asc  —  TC3  0.6  kg  C-4  3.8  m  agl  —  30-152  m 


TX02004.asc  --  TC3  0.6  kg  C-4  3.8  m  agl  —  30-152  m 


TX02007.asc  --  TC2  0.6  kg  C-4  0.3  m  agl  --  91-213  m 


TX02009.asc  --  TC2  0.6  kg  C-4  1.9  m  agl  --  91-213  m 


TX02011.asc  --  TC4  0.6  kg  C-4  1.9  m  agl  --  159-281  m 


Time ,  sec  Time^  sec  Time,  sec 


TX02013.asc  --  TC2  0.6  kg  C-4  1.9  m  agl  —  91-213  m 


Time,  sec  Tinne,  sec  Time,  sec 


TX02015.asc  --  TC4  0.6  kg  C-4  1.9  m  agl  —  159-281  m 


Time,  sec  Time,  sec  Time,  sec 


TX02017.asc  --  TC2  0.6  kg  C-4  1.9  m  agl  --  91-213  m 


Time,  sec  Time,  sec  Time,  sec 


TX02019.asc  --  TC4  0.6  kg  C-4  1.9  m  agl  --  159-281  m 


J\me,  sec  Time,  sec  Time,  sec 


TX02021.asc  --  TC2  2.3  kg  C-4  1.9  m  agl  —  91-213  m 


TX02022.asc  —  TC3  2.3  kg  C-4  1.9  m  agl  ~  30-152  m 


Time,  sec  Time,  sec  Time,  sec 


TX02023.asc  —  TC4  2.3  kg  C-4  1.9  m  agl  —  159-281  m 


Time,  sec  Time,  sec  Time,  sec 


TX02024.asc  —  TCI  2.3  kg  C-4  1.9  m  agl  —  445-567  m 


Time,  sec  Time,  sec  Time,  sec 


TX02025.asc  --  TC2  2.3  kg  C-4  1.9  m  agl  --  91-213  m 


TX02026.asc  --  TC3  2.3  kg  C-4  1.9  m  agl  —  30-152  m 


Q 


Time,  sec  Time,  sec  Time,  sec 


TX02027.ase  --  TC4  2.3  kg  C-4  1.9  m  agl  --  159-281  m 


Time,  sec  Time,  sec  Time,  sec 


TX02029.asc  --  TC2  2.3  kg  C-4  1.9  m  agl  —  91-213  m 


TX02030.asc  --  TC3  2.3  kg  C-4  1.9  m  agl  —  30-152  m 


Time,  sec  Time,  sec  Time,  sec 


TX02031.asc  --  TC4  2.3  kg  C-4  1.9  m  agl  --  159-281  m 


TX02033.asc  --  TC2  2.3  kg  C-4  3.8  m  agl  —  91-213  m 


TX02034.asc  --  TC2  2.3  kg  C-4  3.8  m  agl  --  91-213  m 


Time,  sec  Time,  sec  Time,  sec 


TX02035.asc  --  TC2  2.3  kg  C-4  3.0  m  agl  --  91-213  m 


TX02036.asc  —  TC2  2.3  kg  C-4  1.2  m  agl  --  91-213  m 


Time,  sec  Time,  sec  Time,  sec 


TX02037.asc  --  TC2  2.3  kg  C-4  0.3  m  agl  --  91-213  m 


TX02044.asc  —  TC2  0.6  kg  C-4  3.8  m  agl  —  91-213  m 


TX02045.asc  —  TC2  0.6  kg  C-4  3.0  m  agl  —  91-213  m 


Time,  sec  Time,  sec  Time,  sec 


TX02046.ase  --  TC2  0.6  kg  C-4  1.2  m  agl  —  91-213  m 


Time,  sec  Time,  sec  Time,  sec 


TX02101.ase  —  30m  Pistol  1.0  m--  2-122  m 


TX02102.asc  —  30m  Pistol  1.0  m--  2-122  m 


Time,  sec  Time!  sec  Time,  sec 


TX02103.asc  —  30m  Pistol  1.0  m —  2-122  m 


TX02  —  TX02105.asc  —  90m  Pistol  1.0  m--  2-  62  m 
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Time,  sec  Tima,  sec  Tima,  sec 


TX02  —  TX02106.asc  —  90m  Pistol  1.0  m--  2- 62  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02107.asc  —  90m  Pistol  1.0  m--  2- 62  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02108.asc  —  90m  Pistol  1.0  m--  2- 62  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02109.asc  —  90m  Pistol  1.0  m--  2- 62  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02110.asc  —  60m  Pistol  1.0  m--  2- 92  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02111.asc  --  60m  Pistol  1.0  m--  2- 92  m 
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Tima,  sec  Tima,  sec  Tima,  sec 


TX02  —  TX02112.asc  —  60m  Pistol  1.0  m--  2- 92  m 
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Time,  sec  Time ,  sec  Time,  sec 


Om  Pistol  1.0  m —  2- 92  m 
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Time,  sec  T\me,  sec  Time,  sec 


TX02  —  TX02114.asc  —  60m  Pistol  1.0  m--  2- 92  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02115.asc  —  30m  Pistol  1.0  m—  2-122  m 
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Time,  sec  Time,  sec  Time,  sec 


TX02116.asc  —  30m  Pistol  1.0  m--  2-122  m 


Time,  sec  Time,  sec  Time,  sec 


TX02117.asc  —  30m  Pistol  1.0  m —  2-122  m 


TX02118.asc  —  30m  Pistol  1.0  m —  2-122  m 


Time,  sec  Time,  sec  Time,  sec 


TX02  —  TX02119.asc  —  30m  Pistol  1.0  m--  2-122  m 


Seismic  and  Acoustic  Signature  Measurements  Through  Forest 


89 


Time,  sec  Time,  sec  Time,  sec 
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APPENDIX  B:  MATLAB  PROGRAMS  TO  READ  AND  PLOT 
BINARY  DATA  FROM  NZ  SEISMOGRAPH  AND  TO  WRITE 
MULTICHANNEL  ASCII  DATA  FILES. 

This  section  provides  a  listing  of  the  MATLAB  programs  used  to  read  and 
plot  the  24  binary  data  files  recorded  by  the  Geometries  NZ  seismograph  in  SEG- 
2  format.  The  programs  convert  the  seismograph  voltages  to  physical  units  before 
plotting.  The  programs  also  correct  errors  in  the  shot  time  (due  to  a  malfunction 
of  the  shot  box)  and  write  the  corrected  data  to  multichannel  ASCII  data  files  that 
can  be  plotted  using  the  programs  listed  in  the  next  section  of  this  report.  These 
ASCII  files  serve  as  the  permanent  archive  format  for  the  test  data.  The  files  that 
are  listed  in  this  section  are: 

1 .  tx02doprocessnz.m — Main  program  to  read  binary  data  files  and  write 
multichannel  ASCII  files.  To  use,  set  the  desired  binary  file  numbers  in 
the  variable  “recs”  and  run  the  program.  You  may  have  to  change  the 
lines  that  construct  the  filenames  to  use  the  correct  directory  for  your 
computer,  and  you  may  want  to  change  the  print  options  at  the  very  end 
of  this  program  to  match  your  printer  or  delete  the  plot  from  the  screen 
after  plotting.  The  program  will  write  an  ASCII  file  with  a  new  record 
number  that  either  matches  the  CERL  shot  number  (for  C4  explosions)  or 
a  new  pistol  shot  number  between  101  and  119.  Writes  the  pressure  sen¬ 
sor  info  to  text  file  TX02MikeLog.m.  Calls  the  programs  listed  below. 

2.  readnz.m — Function  to  read  the  SEG-2  binary  file  produced  by  the  NZ 
seismograph.  Calls  words.m.  Writes  info  to  the  text  file  readnzlog.m. 

3.  tx02cerlnum.m — Function  to  select  the  CERL  C4  shot  number  or  new 
pistol  shot  number  from  the  binary  file  number. 

4.  tx021abel.m — Function  to  get  all  of  the  experimental  details  including 
which  sensor  was  attached  to  which  channel,  voltage-to-physical-units 
conversion  factors,  sensor  and  shot  point  geometry  and  distances,  which 
source  used  for  each  file,  etc.  Calls  words.m. 

5.  tx02plotnz.m — Function  to  make  a  three-panel  plot  of  the  binary  data  in 
landscape  orientation.  The  channels  plotted  in  each  of  the  panels  are 
given  by  the  contents  of  iplotl,  iplot2,  and  iplot3.  These  plots  show  the 
uncorrected  times  in  the  binary  files,  i.e.,  when  the  recorder  button  was 
pushed,  not  when  the  shot  was  actually  fired. 

6.  tx02tbcorr.m — Function  to  find  the  true  time  break  (shot  instant).  This 
is  done  using  the  arrival  times  from  the  first  and  last  sensor  in  the  array 
to  determine  the  acoustic  wave  velocity,  and  then  projecting  back  to  zero 
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distance  to  get  the  origin  time.  Writes  the  time  break  info  to  file 
Tx02TimeBreakLog.m  and  applies  the  correction  to  the  multichannel 
ASCII  data  files. 

7.  tx02writenchan.m — Function  to  write  multichannel  ASCII  data  file 
with  12  line  trace  headers.  The  ASCII  file  name  is  tx02nnn.asc,  where 
nnn  is  a  new  file  number,  nnn  =  CERL  shot  number  if  it  was  a  C4  shot. 

8.  words.m — Function  to  separate  a  string  into  individual  words  by  finding 
the  blank  characters. 

tx02doprocessnz.m — Main  program  to  plot  binary  data  and  write  multichannel  ASCII  file 

%  d: \DataWinNT\FY02LoneStarTX\Plot\tx02doprocessnz .m 
%  Processes  TX02  NZ  seismograph  data: 

%  Reads  NZ  data  (x)  --  writes  file  info  to  readnzlog.m 

%  Converts  to  eng  units  (xx) 

%  Finds  CERL  rec  number  (nnn) 

%  Plots  data  -  3  panels 

%  Finds  the  true  time  break  -  writes  in  Tx02TimeBreakLog . m 
%  Writes  pressure  data  info  to  Tx02MikeLog . m 
%  Writes  multichannel  ascii  file  to  tx02nnn.asc 
%  Calls  subroutines 
%  readnz  words  tx02cerlnum 
%  tx02plotnz  tx021abel  tx02tbcorr 
%  tx02writenchan 

%  d:\DataWinNT\FY02LoneStarTX\Plot\tx02doprocessnz.m 
%  Reads  NZ  seismograph  binary  records,  makes  3  panel  plots  of  data 
%  Writes  pressure  sensor  info  to  Tx02MikeLog . m 

o 

o 

%  Calls  functions 

%  readnz. m  -  to  read  the  SEG-2  binary  data  using  the 
%  CRREL  field  record  numbers 

%  words.m  -  to  parse  ascii  header  lines  into  words 
%  Writes  info  about  the  file  to  file  readnzlog.m 

%  tx02cerlnum . m  -  selects  CERL  shot  number  or  new  pistol  shot  number 
%  calls  no  subroutines 

%  tx021abel.m  -  to  get  sensor  geometry,  etc 
%  calls  words.m 

%  tx02plotnz.m  -  3  Panel  plot  of  data 
%  calls  no  subroutines 

%  tx02tbcorr.m  -  Finds  true  time  break  (shot  instant) 
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%  from  the  arrival  times  at  a  few  of  the  sensors 
%  Writes  info  to  Tx02TimeBreakLog . m 
%  calls  no  subroutines 

%  tx02writenchan . m  -  writes  multichannel  ascii  data  with  12  line  trace 
headers 

%  File  name  is  tx02nnn.asc,  where  nnn  is  a  new  file  number 
%  nnn  =  CERL  Shot  number  if  it  was  a  C-4  shot 
%  calls  no  subroutines 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 


%  NZ  binary  record  numbers;  plot  has  new  record  number  in  title 
%  C-4  Records  39,44  were  bad  and  are  omitted 

%recs  =  [65  66  67  68  5];  %TX02  C-4  examples  and  blank  pistol  shot 
example 

%recs  =  [  3  5  10  15] ;  %Blank  pistol  shot  examples 

%recs  =  [  40:43  60:68  70:74  93:113  1:19];  %A11  C-4  39,44  bad  and  all 
blank  pistol 

recs  =  [  40:43  60:68  70:74  93:113  1:19];  %A11  C-4  39,44  bad  and  all 
blank  pistol 
recs  =  [64  65 ]  ; 

%  Loop  to  plot  the  files 

for  i  =  1 : length (recs) ;  %  **************************************** 
xx=[]  ;x=[]  ; 
recnum  =  recs  (i) ; 
shotnum  =  recs  (i) ; 

fname  =  sprint f (  '  d: \\DataWinNT\\f y02texas\\%s . dat 1 , num2str (recnum) ) ; 

%  read  in  data  from  NZ  binary  file  -  x 
[x, npts , deltat , nchan, delayms , descalingf act ,  . . . 
stackcount , acqdate, acqtime]  =  readnz (fname) ; 

%  sample  rate 

srate  =  1/  deltat ; %number  of  samples  per  second 

%  Get  geometry  and  description  for  this  experiment  and  record  number 
[sensor, sensortype, rtype, sensorx, sensory, sensorz, . . . 
srcloc, sourcetype, stype, srcx, srcy, srcz, sourcesize, . . . 
engf act , dist , sensor2 , source2 ]  =  tx021abel  (x, shotnum, nchan) ; 
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%  Convert  from  voltage  to  physical  units  (Pa  or  m/s) 
xgain  =  descalingf act . / (stackcount . *engf act ) ;  %recorder  units  to  eng 
units 

for  iii  =  1  inchan;  xx(:,iii)  =  x(:,iii)  *  xgain (iii);  end 
xxmax  =  max (abs (xx) ) ;  %abs  max  for  plotting 

%  get  CERL  record  number  for  C-4  shots 
cerlrec  =  0; 

[cerlrec  ]  =  tx02cerlnum  (shotnum) ; %If  C-4  get  CERL  Shot  Number  for 
Plot 

%  3  panel  plot  of  nz  data 

iplotl  =  [ 9  10  11  12  5  6  7  8  1  2  3  4 ] ;  %Channels  for  1st  (left)  panel 
iplot2  =  [22  21  23  24  18  17  19  20  14  13  15  16  ];%  2nd  (center)  panel 
iplot3  =  [24  20  16  11  12  7  8  3  4];%Channels  for  3rd  (right)  panel 
tx02plotnz ( shotnum, cerlrec , xx , iplotl , iplot2 , iplot3 ,  del tat ,  delayms , . . . 
dist,  sensor,  sensorx, sensory, sensorz,  .  .  . 
srcloc, sourcesize, sourcetype, srcx, srcy, srcz) ; 

%  Write  pressure  sensor  data  to  text  file 
f id3=f open ( 1 Tx02MikeLog . m 1 ,  '  a  1 ) ; 
iplot  =  iplot3;  nplot=length ( iplot3 ) ; 

%  Write  pressure  sensor  trace  header  info  to  logfile 
[aa  bb]  =  max (xx) ;  %set  plot  windows  because  manual  TB 
[cc  dd]  =  min (xx) ;  %set  plot  windows  because  manual  TB 
for  iii  =l:nplot; 

tttmax  =  delayms/1000  +  (bb  ( iplot  ( iii ) )  ) /srate;  %times  in  seconds 

tttmin  =  delayms/1000  +  (dd ( iplot ( iii ) )  ) /srate;  %times  in  seconds 

fprintf (fid3, . . . 

'  %4 . Of  %3 . Of  %3 . Of  %3 . If  %5 . Of  %3.1f  %3.1f  %6.1f  %6.1f  %6.4f  %6.1f 
%6 . 4 f  %6.4f  \n' ,  .  .  . 

shotnum, iplot (iii) , sensorx (iplot (iii) ) , sensorz (iplot (iii) ) , . . . 
srcx, srcz , sourcesize*l . 25/2 . 2 ,  .  .  . 

dist (iplot (iii) ) , aa (iplot (iii) ) , tttmax, cc (iplot (iii) ) , tttmin, (tttmin- 
tttmax) ) ; 
end 

f close ( f id3 ) ; 


o 

o 


Find  TB  correction  -  writes  to  output  file  xxxTimebreaklog . m 
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%  itimeTB  is  the  TB  (shot  instant)  index 
%  itimeO  is  the  index  to  start  the  save  (usually=TB) 

%  if  the  shot  was  at  TCI,  1  second  is  added  to  itimeO 
%  itimel  is  the  index  2  seconds  after  the  TB 
%  timeO  is  the  start  time  for  saving  the  time  series 

[ itimeTB, itimeO , itimel , timeO , xdist, xtime]  =  tx02 tbcorr ( shotnum, . . . 
cerlrec, xx, del tat, delayms , srcloc) ; 

%  Write  engineering  data  to  ascii  file 

tx02writenchan (shotnum, cerlrec, xx, npts, deltat, nchan, delayms, descalingfac 

t,  .  .  . 

stackcount, acqdate, acqtime, engf act , xgain, itimeO , itimel , timeO , . . . 
dist,  sensor,  rtype, sensorx, sensory, sensorz, sensor2,  . . . 
srcloc, stype, sourcetype, sourcesize, srcx, srcy, srcz,  source2) ; 

%  Print  and  delete 
%  print  -dljetplus; 

%  print  -dljetplus;  pause (25);  h=gcf;  delete (h) ; 

%  print  to  PostScript  file  for  conversion  to  PDF: 

%  print  -dps  -rl50  -append  tx02datafigs 
%  h=gcf; delete (h) ; 

end  %loop  over  files  %  **************************************** 
f close ( 'all1 ) ; 

return 

tx02label.m — Provides  info  about  sensor  array  and  sources 

function  [sensor, sensortype, rtype, sensorx, sensory, sensorz, . . . 
srcloc, sourcetype, stype, srcx, srcy, srcz, sourcesize, . . . 
engf act , dist , sensor2 , source2 ]  =  tx021abel  (x, shotnum, nchan) 

%  d:\DataWinNT\FY02LoneStarTX\tx021abel.m 

%  Returns  channel  ids  for  NZ  seismograph  used  at  FY02  TX  test 
%  calls  words. m  -  breaks  string  into  individual  words 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 
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%  Assign  sensor  coordinates 

%  sensorx  =  ( 1 : nchan) *1 ;  sensorz ( 1 : nchan)  =  0; 

sensorx  =  (1 :nchan) *0;  sensory  =  (1 :nchan) *0;  sensorz  =  (1 :nchan) *0; 
sensorx  =  [ 30 , 30 , 30 , 30 ,  60,60,60,60,  90,90,90,90,  ... 

120,120,120,120,  150,150,150,150,  152,152,152,152] ; 

heightl5  =  [4,8,12,24];  sensorz (heightl5)  =  1.5; 
height2  =  [3,7,11];  sensor z (height2 )  =  0.08; 
height3  =  [16,20];  sensorz (height3)  =  0.22; 

%  Assign  sensortype 

for  i  =  1  inchan;  sensortype ( i )  =  ' P 1 ;  end 
senv  =  [2  6  10  15  19  23  ];  sensortype ( senv) =' V 1 ; 
senr  =  [1  5  9  13  17  21  ];  sensortype ( senr) =' R ' ; 
sent  =  [  14  18  22  ];  sensortype (sent )=' T 1 ; 

sss  =  [ ' R-30  V-30  P-30-0  P-30-1  R-60  V-60  P-60-0  P-60-1  '  ... 

' R- 90  V-90  P-90-0  P-90-1  '  ... 

' R-12 0  T-120  V-120  P-120-0  R-150  T-150  V-150  P-150-0  '  ... 

' R-152T  T-152T  V-152T  P-152T-1  ' ]; 

% s s s  =  [ ' 1  2  3  4  5  6  7  8  9  10  11  12  ’  ... 

%  ’13  14  15  16  17  18  19  20  21  22  23  24  1  ... 

%  '25  26  27  28  29  30  31  32  33  34  35  36  '  ]  ; 

sensor  =  words (sss)  ; 

sss2  =  [ 'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-102A07-15972  PCB- 
102A07-15973  '  ... 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-10 6B50-32 60  PCB-10 6B50-6693  ' 

'MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  PCB-106B50-3259  PCB-10 6B50-6695  ' 

1 MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4939-1  '  .. 

' MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  MP-4 . 5Hz-geophone  B&K-4939-2  '  .. 

' MP-4 . 5Hz-geophone-in-tree  MP-4 . 5Hz-geophone-in-tree  '  ... 

1 MP-4 . 5Hz-geophone-in-tree  B&K-4 93 9-3-in-tree  ']; 
sensor2  =  words (sss2); 

%  rtype  1  =  Microphone,  2  =  PCB,  5=  Vert,  6  =  Horiz,  7  =  Accel 
rtype  =  [652265226522665166516651]; 
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%  engfact  =  units  per  mV 

%  mV  =  data_value*descal_f act/stack_count 

%  phys_units  =  data_value*descale_f act  /  (stack_count*engf act ) 


%  Starting  nominal  values : 


engfact  =  [  . . . 


32.2*1000 

32.2*1000 

32.2*1000 

32.2*1000 

32.2*1000 

32.2*1000 


32.2*1000 

32.2*1000 

32.2*1000 

32.2*1000 

32.2*1000 

32.2*1000 


100/6895  100/6895  .  .  . 
500/6895  500/6895  .  .  . 
500/6895  500/6895  .  .  . 
32.2*1000  -1.2/1000  .  . 
32.2*1000  -1.2/1000  .  . 
32.2*1000  -1.2/1000  ]  ; 


if (shotnum  >=45) ; 

engfact  (7)  =  100/6895;  sensor2  (7 ,  : )  =  1 PCB102A07-15971  ' 
end  %102A07  15971  from  Rec  45  on 


%  Actual  calibrated  values: 
engfact  =  [  ... 

32.2*1000  32.2*1000  88.09/6895  85.75/6895  ... 

32.2*1000  32.2*1000  500/6895  464.5/6895  ... 

32.2*1000  32.2*1000  563.2/6895  500/6895  ... 

32.2*1000  32.2*1000  32.2*1000  -1.15/1000  ... 

32.2*1000  32.2*1000  32.2*1000  -1.32/1000  ... 

32.2*1000  32.2*1000  32.2*1000  -1.20/1000  ]; 

%  Adjust  for  changes  in  sensor  array  during  test 

if (shotnum  >=45);  engfact (7)  =  89.72/6895;  end  %102A07  15971  from  Rec  45 
on 

if (shotnum  >=5  &  shotnum  <=27);  %B&K  mikes  used  for  blank  pistol  shots 
engfact (16)  =-0.7501/1000;  sensor2 ( 1 6 , : ) = 1 B&K-4 1 65-5  '; 
engfact(20)  =  -0.7436/1000;  sensor2  (20,  : )  =  ' B&K-4165-6  '; 
end  % 

%  Assign  source  properties 
%  Default  values 

srcloc  =  1 SP1 1 ;  %Nominal  value,  actual  TC  locations  below 
srclocnum  =  1;  srcx  =  0.0;  srcy  =  0.0;  srcz  =  1.5; 

%  stype  1  =  C-4,  6  =  50  cal,  9  =  45  cal,  10  = 
sourcetype  =  'C-4’;  stype  =  1;  sourcesize  =  % ; 


noise 
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%  List  of  source  heights  used  during  tests 

shtl  =  0.3048*  (6*12  +  3) / 1 2 ; 

sht2  =  0.3048* (12*12+6) /12; 

sht3  =  0.3048* (10*12+0) /12; 

sht4  =  0.3048*  (4*12  +  0) /12; 

sht5  =  0.3048*  (2*12  +  0) /12; 

sht5  =  0.3048* (1*12+0) /12; 

srcz  =  shtl;  %default  source  height  =6  ft  3=  1.9m 
%  Assign  C-4  source  size 

if (shotnum==40  |  (shotnum>=96  &  shotnum<=113 )  )  %4  bricks 

sourcesize  =  4; 
end 

%  Assign  source  height 

if (shotnum==41  |  shotnum==70  |  shotnum==109  ) 
srcz  =  sht2;  %12-6  feet  and  inches 
end 

if (shotnum==42  |  shotnum==71  |  shotnum==110  ) 
srcz  =  sht2;  %10-0 
end 

if (shotnum==43  |  shotnum==72  |  shotnum==lll  ) 
srcz  =  sht3;  %4-0 
end 

if (shotnum==44  |  shotnum==73  |  shotnum==112  ) 
srcz  =  sht4;  %2-0 
end 

if (shotnum==60  |  shotnum==74  |  shotnum==113  ) 
srcz  =  sht5;  %l-0 
end 


%  Assign  source  locations 

%  TC3  is  the  origin,  30  m  from  the  first  CRREL  sensor 
if (shotnum==61  |  shotnum==65  |  shotnum==65  . . . 

|  shotnum==93  |  shotnum==96  |  shotnum==100  . . . 

|  shotnum==104  )  %TC1 

srclocnum=l ; srcloc= ' TCI ' ; srcx=-414 . 7 ; source 2= ' Open-field ' ; 
end 

if (shotnum==60  |  shotnum==62  |  shotnum==66  . . . 

|  (shotnum>=70  &  shotnum<=74)  |  shotnum==94  ... 
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|  shotnum==97  |  shotnum==101  |  shotnum==105  . . . 

|  (shotnum>=108  &  shotnum<=113 )  )  %TC2 

srclocnum=2 ; srcloc= ' TC2 ' ; srcx=-60 . 5 ; source 2= 1 0pen-field ' ; 
end 

if ( (shotnum>=39  &  shotnum<=44)  |  shotnum==63  . . . 

|  shotnum==67  |  shotnum==98  |  shotnum==102  . . . 

|  shotnum==106  )  %TC3 

srclocnum=3 ; srcloc= ' TC3 ' ; srcx=0 . 0 ; source 2= ' Edge-of- forest ' ; 
end 

if (shotnum==64  |  shotnum==68  . . . 

|  shotnum==95  |  shotnum==99  |  shotnum==103  . . . 

|  shotnum==107  )  %TC4 

srclocnum=4 ; srcloc= 1 TC4 ' ; srcx=311 . 0 ; source 2= ' Edge-of- forest ' ; 
end 

if (shotnum==39  |  shotnum==69) 
srclocnum=5 ; srcloc= '  ' ; srcx=0 . 0 ; 

sourcetype  =  'Noise’;  stype  =  10;  sourcesize  =  0 ; source2= ' Noise ' ; 
end 

%  Assign  source  locations  for  blank  pistol  shots 
if (shotnum<5)  %30m 

srclocnum=5 ; srcloc= ' 30m ' ; srcx=30 . 0 ; srcy=2 ;srcz=l; 

sourcetype  =  'Pistol';  stype  =  9;  sourcesize  =  0 ; source2= ' In-forest ' ; 
end 

if (shotnum>=5  &  shotnumClO)  %90m 
srclocnum=7 ; srcloc= ' 90m ' ; srcx=90 . 0 ; srcy=2 ; srcz=l; 

sourcetype  =  'Pistol';  stype  =  9;  sourcesize  =  0 ; source2= ' In-forest ' ; 
end 

if (shotnum>=10  &  shotnum<15)  %60m 
srclocnum=6 ; srcloc= ' 60m ' ; srcx=60 . 0 ; srcy=2 ;srcz=l; 

sourcetype  =  'Pistol';  stype  =  9;  sourcesize  =  0 ; source2= ' In-forest ' ; 
end 

if (shotnum>=15  &  shotnum<20)  %30m 
srclocnum=5 ; srcloc= ' 30m ' ; srcx=30 . 0 ; srcy=2 ; srcz=l; 

sourcetype  =  'Pistol';  stype  =  9;  sourcesize  =  0 ; source2= ' In-forest ' ; 
end 


%  Calculate  source-to-sensor  distances  for  this  shot 
srcx2  =  srcx*srcx;  srcy2  =  srcy^srcy;  srcz2  =  srcz^srcz; 
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for  i=l :nchan 

dist(i)  =  ( (sensorx (i) -srcx) . A2  +  ... 

(sensory (i) -srcy) . A2  +  ... 

( sensorz  ( i ) -srcz )  . A2  )  .A0.5; 

end 

return  % - 

tx02cerlnum.m — Changes  CRREL  binary  file  numbers  to  CERL  shot  numbers 

function  [cerlrec  ]  =  tx02cerlnum  (shotnum) 

%  d:\DataWinNT\FY02LoneStarTX\tx02cerlnum.m 

%  Converts  CRREL  NZ  seismograph  file  numbers  for  FY02  TX  test 
%  For  C-4,  returns  CERL  Shot  Num  (1-39,  with  some  "B"  repeats) 
%  For  blank  pistol,  assigns  a  number  in  the  100's 
%  For  large  Demo  shots,  assigns  a  number  in  the  300's 
%  Noise  records  92->  401 

%  Cal  record  numbers  are  omitted  20-27,  82-91,  114-199 
%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 
cerlrec  =  '  ' ; 

if ( shotnum==3 9 ); cerlrec  =  '  1 ' ; end 

if ( shotnum==4 0 ); cerlrec  =  '  2 ' ; end 

if ( shotnum==4 1 ); cerlrec  =  '  3 ' ; end 

if ( shotnum==42 ); cerlrec  =  '  4 ' ; end 

if (shotnum==43 ); cerlrec  =  '  5 ' ; end 

if ( shotnum==4 4 ); cerlrec  =  '  6 ' ; end 

if ( shotnum==60 ); cerlrec  =  '  7 ' ; end 

if ( shotnum==61 ); cerlrec  =  '  8 ' ; end 

if ( shotnum==62 ); cerlrec  =  '  9 ' ; end 

if ( shotnum==63 ); cerlrec  =  '  10'; end 

if ( shotnum==64 ); cerlrec  =  '  11'; end 

if ( shotnum==65 ); cerlrec  =  '  12'; end 

if (shotnum==66) ; cerlrec  =  '  13'; end 

if ( shotnum==67 ); cerlrec  =  '  14'; end 

if ( shotnum==68 ); cerlrec  =  '  15'; end 

%  The  next  shots  are  repeats  of  Shots  3-7 

%  We  assign  the  CERL  numbers  43-47 

if ( shotnum==7 0 ); cerlrec  =  '  43'; end 

if ( shotnum==7 1 ); cerlrec  =  '  44'; end 
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if ( shotnum==72 ) ; cerlrec  =  '  45'; end 
if ( shotnum==73 ) ; cerlrec  =  '  46'; end 
if ( shotnum==7 4 ); cerlrec  =  '  47'; end 
if ( shotnum==93 ); cerlrec  =  '  16'; end 
if ( shotnum==94 ); cerlrec  =  '  17'; end 
if ( shotnum==95 ); cerlrec  =  '  19'; end 
%  missed  CERL  shot  18 
if (shotnum==96) ; cerlrec  =  '  20'; end 
if ( shotnum==97 ); cerlrec  =  '  21'; end 
if ( shotnum==98 ); cerlrec  =  '  22'; end 
if (shotnum==99) ; cerlrec  =  '  23'; end 
if (shotnum==100 ); cerlrec  =  '  24'; end 
if (shotnum==101 ); cerlrec  =  '  25'; end 
if (shotnum==102 )/ cerlrec  =  '  26'; end 
if (shotnum==103) / cerlrec  =  '  27'; end 
if (shotnum==104 ); cerlrec  =  '  28'; end 
if (shotnum==105 )/ cerlrec  =  '  29'; end 
if (shotnum==106) ; cerlrec  =  ’  30'; end 
if (shotnum==107 ); cerlrec  =  '  31'; end 
if (shotnum==108 ); cerlrec  =  '  32'; end 
if (shotnum==109) ; cerlrec  =  '  33'; end 
if (shotnum==110 ); cerlrec  =  '  34'; end 
if (shotnum==lll ); cerlrec  =  '  35'; end 
if ( shotnum==l 12 ); cerlrec  =  '  36'; end 
if (shotnum==113 ); cerlrec  =  ’  37'; end 
%  CRREL  Blank  Pistol  Shots 
if (shotnum==  1); cerlrec  =  '101’; end 
if (shotnum==  2); cerlrec  =  '102'; end 
if (shotnum==  3) /cerlrec  =  '103'; end 
if (shotnum==  4) /cerlrec  =  '104' /end 
if (shotnum==  5) /cerlrec  =  '105' /end 
if (shotnum==  6) /cerlrec  =  '106' /end 
if (shotnum==  7) /cerlrec  =  '107' /end 
if (shotnum==  8) /cerlrec  =  '108' /end 
if (shotnum==  9) /cerlrec  =  '109' /end 
if (shotnum==10 )/ cerlrec  =  '110' /end 
if (shotnum==ll )/ cerlrec  =  '111' /end 
if ( shotnum==12 )/ cerlrec  =  '112' /end 
if ( shotnum==13 )/ cerlrec  =  '113' /end 
if (shotnum==14 )/ cerlrec  =  ’114’ /end 
if ( shotnum==15 )/ cerlrec  =  '115' /end 
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if ( shotnum==l 6 ) ; cerlrec  =  '116'; end 
if ( shotnum==17 ) ; cerlrec  =  '117'; end 
if (shotnum==18 ); cerlrec  =  '118'; end 
if ( shotnum==l 9 ); cerlrec  =  '119'; end 
%  Large  in-ground  Demolition  Shots 
if ( shotnum==2 8 ); cerlrec  =  '301'; end 
if ( shotnum==2 9 ); cerlrec  =  '302'; end 
if ( shotnum==30 ); cerlrec  =  '303'; end 
if ( shotnum==31 ); cerlrec  =  '304'; end 
if ( shotnum==32 ); cerlrec  =  '305'; end 
if ( shotnum==33 ); cerlrec  =  '306'; end 
if ( shotnum==34 ); cerlrec  =  '307'; end 
if ( shotnum==35 ); cerlrec  =  '308'; end 
if ( shotnum==3 6 ); cerlrec  =  '309'; end 
if ( shotnum==37 ); cerlrec  =  '310'; end 
if ( shotnum==38 ); cerlrec  =  '311'; end 
if (shotnum==45 ); cerlrec  =  '312'; end 
if ( shotnum==4 6 ); cerlrec  =  '313'; end 
if (shotnum==47 ); cerlrec  =  '314’; end 
if ( shotnum==4 8 ); cerlrec  =  '315'; end 
if ( shotnum==4 9 ); cerlrec  =  '316'; end 
if ( shotnum==50 ); cerlrec  =  '317'; end 
if ( shotnum==51 ); cerlrec  =  '318'; end 
if ( shotnum==52 ); cerlrec  =  '319'; end 
if ( shotnum==53 ); cerlrec  =  '320'; end 
if ( shotnum==54 ); cerlrec  =  '321'; end 
if ( shotnum==55 ); cerlrec  =  '322'; end 
if ( shotnum==5 6 ); cerlrec  =  '323'; end 
if ( shotnum==57 ); cerlrec  =  '324'; end 
if ( shotnum==58 ); cerlrec  =  '325'; end 
if ( shotnum==5 9 ); cerlrec  =  '326'; end 
if ( shotnum==75 ); cerlrec  =  '327'; end 
if ( shotnum==7 6 ); cerlrec  =  '328'; end 
if (shotnum==77 ); cerlrec  =  '329'; end 
if ( shotnum==7 8 ); cerlrec  =  '330'; end 
if ( shotnum==7 9 ); cerlrec  =  '331’; end 
if (shotnum==80 ); cerlrec  =  '332'; end 
if (shotnum==81 ); cerlrec  =  '333'; end 
%  Misfire  =  Noise  record 
if ( shotnum==92 ); cerlrec  =  '401'; end 
return 
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tx02plotnz.m — Three-panel  plot  of  binary  seismograph  data 

function 

tx02plotnz (shotnum, cerlrec, xx, iplotl, iplot2, iplot3,  deltat,  delayms,  . . . 
dist,  sensor ,  sensorx, sensory, sensorz ,  . . . 
srcloc, sourcesize, sourcetype, srcx, srcy, srcz) ; 

%  d:\DataWinNT\FY02LoneStarTX\Plot\tx02plotnz.m 

o 

o 

%  Plots  nchan  channels  of  data  for  Texas  02  NZ  data 
%  Data  are  in  array  xx 

%  Channels  to  plot  in  each  panel  are  in  iplotl, 2, 3 
%  Does  not  call  any  subroutines 

o 

o 

%  Because  of  CRREL  blaster  box  problems,  there  are  no  good  time 
%  breaks  in  the  data,  and  long  record  lengths  were  used.  The 
%  data  are  aligned  before  plotting  based  on  the  actual  arrival 
%  times  found  in  the  recorded  data. 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

nplotl  =  length ( iplotl ); nplot2  =  length ( iplot2 ) ; 
nplot3  =  length ( iplot3 ) ; 

xxmax  =  max (abs (xx) ) ;  %abs  max  for  plotting 
srate  =  1/  deltat;  %number  of  samples  per  second 


%  Determine  the  time  window  and  direction  of  arrivals  for  the  shots 
%  Time  windows  for  C-4  shots 

[aa  bb]  =  max (xx) ;  %set  plot  windows  because  manual  TB 
[cc  dd]  =  min (xx) ;  %for  output  file 

tmax4  =  delayms/1000  +  (bb(4)  ) /srate;  %times  in  seconds 

tmax20  =  delayms/1000  +  (bb  ( 2 0 ) ) /srate ; 

tmaxl6  =  delayms/1000  +  (bb ( 16) ) /srate; 

tmaxl2  =  delayms/1000  +  (bb ( 12 )) /srate; 

tstart4=  floor (tmax4  *10) /10; 

tstart20=  round (tmax20*10 ) /10 ; 

tstartl6=  f loor ( tmaxl 6*10 ) /10 ; 

if  (tstart4  <  tstartl6) 
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tstartpanell  =  tstart4 ; tendpanell  =  tstartpanell  +  0.5; 
tstartpanel2  =  tstartl 6 ; tendpanel2  =  tstartpanel2  +  0.5; 
tstartpanel3  =  tstartpanell ; tendpanel3  =  tstart20  +  0.1; 
else  %if (tstartl6  <  tstart4) 
tstart4=  round(tmax4  *10) /10; 
tstartl2=  f loor (tmaxl2*10 ) /10 ; 
tstart20=  f loor (tmax20*10 ) /10 ; 
tstartl6=  round (tmaxl 6*10 ) /10 ; 
tstartpanell  =  tstartl2 ; tendpanell  =  tstartpanell  +  0.5; 
tstartpanel2  =  tstart2 0 ; tendpanel2  =  tstartpanel2  +  0.5; 
tstartpanel3  =  tstartpanel2 ; tendpanel3  =  tstart4  +  0.1; 
end 

%  Time  windows  for  blank  pistol  shots 

if (shotnum<5)  %Pistol  shots  at  30  m,  use  mike  at  90  m  =  Chan  12 
tstartl2=  floor (tmaxl2*10 ) /10 ; 

tstartpanell  =  tstartl2-0 . 2 ; tendpanell  =  tstartpanell  +  0.6; 

tstartpanel2  =  tstartl2-0 . 2 ; tendpanel2  =  tstartpanel2  +  0.6; 

tstartpanel3  =  tstartpanell ; tendpanel3  =  tendpanell; 

end 

if (shotnum>4  &  shotnum<20)  %Pistol  shots  at  90  m,  use  mike  at  90  m  = 
Chan  12 

tstartl2=  floor (tmaxl2*10 ) /10 ; 

tstartpanell  =  tstartl2-0 . 1 ; tendpanell  =  tstartpanell  +  0.6; 

tstartpanel2  =  tstartl2-0 . 1 ; tendpanel2  =  tstartpanel2  +  0.6; 

tstartpanel3  =  tstartpanell ; tendpanel3  =  tendpanell; 

end 

%  time  series 

tend  =  length (xx (:, 1 )) ;  %all  data 
t  =  delayms/1000  +  ( 1 : tend) /srate ; 

%  Start  of  data  plot  %*************************************** 

nametext  =  [fNZ  Data  -  3  panel'];  % - FIGURE  -  3  Panels 

fl  =  f igure (' Name ' , nametext ) ; 

set (gcf , ' Units ' , ' inches ' , ' PaperUnits ' , ' inches ' ) ; 

subplot ( 1 , 3 , 1 ) ;  %Left  (1st)  panel  of  plot  - 

iplot=iplotl ;  nplot=nplotl ; 

%iplot=l:12;  %Channels  1-12  %iplot=12 : -1 : 1 ;  %Channels  1-12 
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%  Set  axis  params  here;  depends  on  nplot  -  t2 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set  (axl ,  1 YTick 1 ,  [ ] ) ;  set  (axl ,  ' Box  1 ,  ' on ' ) ; 
xlabel ( ' Time,  sec') 

%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

plot  (t,  (xx (1 : length (t) , iplot  (i) ) /xxmax (iplot (i) ) )  +y shift (i) ,  ' kf ) ; 
leftlabel  =  sprintf (' %g round (dist (iplot (i) ))) ; 

%  Text  labels  for  plot 
strl  =  sensor ( iplot ( i ) ,  :) ; 

str2  =  sprintf ( ' %7 . 2e 1 ,  (xxmax (iplot (i) ) ) )  ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( 1 %4 . lfE-6 1 ,  (IE 6* xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 
str2=sprintf ( ' %4 . 2  fE-3 ' ,  ( lE3*xxmax ( iplot ( i ) ) ) )  ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( 1 %4 . lfE-3 ' , (IE 3* xxmax (iplot (i) ) ) ) ; 

elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( ' %4 . 2 f 1 ,  (xxmax ( iplot ( i )  )  )  )  ; 

elseif  (xxmax ( iplot ( i ))  )  >  1000.0 

str2=sprintf (  '  %4 . If k ' , (  xxmax (iplot (i))/1000  )); 

else 

str2=sprintf ( ' %4 . If  1 ,  (xxmax (iplot (i) ) ) ) ; 

end  %text  loop 

str3  =  char (strl); 

str4  =  char(str3  ,  str2); 

str5  =  cellstr (str4 ) ; 

text ( tendpanell  +  0 . 02 ,  [yshift (i) ] , str5,  ' Fonts ize 1 A  8,  .  .  . 

' HorizontalAlignment ' ,  1  left ' ) ; 
end  %plot  loop  i=l : nplot 

subplot ( 1 , 3 , 2 ) ;  %Center  (2nd)  panel  of  plot  - 

iplot=iplot2 ;  nplot=nplot2 ; 
yshift  =  ( 1 : nplot ) *2  -  2; 

ax  =  [ tstartpanel2-0 . 05  tendpanel2  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set (axl ,  'YTick ’ ,  [ ] ) ;  set (axl ,  ' Box ' ,  1  on ' ) ; 
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xlabel ( 1  Time,  sec') 

%  plot  data  and  label 

hold  on; 

for  i=l:nplot 

plot  (t,  (xx (1 : length (t) , ip lot (i) ) /xxmax (iplot (i) ) )  +y shift (i)  ,  ' k’  )  ; 
leftlabel  =  sprintf (' %g round (dist (iplot (i) ))) ; 

%  Text  labels  for  plot 
strl  =  sensor ( iplot ( i ),:) ; 

str2  =  sprintf ( 1 %7 . 2e 1 ,  (xxmax ( iplot ( i ) ) ) ) ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( ' %4 . lfE-6 ' , (lE6*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 
str2=sprintf ( 1 %4 . 2  fE-3 ' ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf (  '  %4 . lfE-3 ' ,  (lE3*xxmax (iplot (i) ) ) )  ; 

elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( 1 %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 

elseif  (xxmax ( iplot ( i )) )  >  1000.0 

str2=sprintf ( ' %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 

else 

str2=sprintf ( 1 %4 . If ' ,  (xxmax (iplot (i) ) ) ) ; 

end  %text  loop 

str3  =  char (strl); 

str4  =  char(str3  ,  str2); 

str5  =  cellstr ( str4  )  ; 

text ( tendpanel2  +  0 . 02,  [y shift (i) ] Astr5,  ' Font Size ',8,  .  .  . 

1 HorizontalAlignment ' , ' left ' ) ; 
end  %plot  loop  i=l:nplot 


subplot ( 1 , 3 , 3 ) ;  %Right  (3rd)  panel  of  plot  - 

%  Pressure  sensors  only 
iplot=iplot3 ;  nplot=nplot3 ; 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanel3-0 . 05  tendpanel3  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set  (axl ,  1 YTick 1 ,  [ ] ) ;  set (axl ,  ' Box ' ,  ' on ' ) ; 
xlabel (' Timef  sec') 

o 
o 
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hold  on; 
for  i=l : nplot 

if (xxmax ( iplot ( i ) ==0 ) ) ; xxmax ( iplot ( i ) ) =1 ; end 

plot  (t,  (xx (1 : length (t) , iplot  (i) ) /xxmax (iplot (i) ) )  +y shift (i)  ,  ’ k 1  ) 
leftlabel  =  sprintf (' %g ' ,  round (dist (iplot (i) )))  ; 

%  Text  labels  for  plot 
strl  =  sensor (iplot (i) ,  :) ; 

str2  =  sprintf ( 1 %7 . 2e 1 , (xxmax (iplot (i) ) ) ) ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( 1 %4 . lfE-6 1 ,  (IE 6* xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 
str2=sprintf ( ' %4 . 2  fE-3 ' , ( lE3*xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( 1 %4 . lfE-3 ' , (IE 3* xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  10.0 
str2=sprintf ( 1 %4 . 2  f 1 ,  (xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  >  1000.0 
str2=sprintf ( ' %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 
else 

str2=sprintf ( ' %4 . If  1 ,  (xxmax (iplot (i) ) ) ) ; 
end  %text  loop 
str3  =  char (strl); 
str4  =  char(str3  A  str2); 
str5  =  cellstr (str4 ) ; 

text ( tendpanel3  +  0 .02,  [y shift (i) ]  ,  str5,  ' Font Size  1 , 8,  .  .  . 

' HorizontalAlignment ' ,  1  left ' ) ; 

%  Write  title  for  right  panel 
if (i==nplot) 

text ( tstartpanel3  +  0 . 05, yshift (i) +1,  ' Pressure  Sensors  1 ) ; 
end 


end  %plot  loop  1=1 : nplot 


%  Title  for  entire  plot 


if ( sourcesize~=0 )  %C-4  shot 
plottitle  =  ... 

sprintf (' TX02  --  Shot  %s  Rec  %g  --  %s  %3.1f  kg  %s  %3.1f  m  agl  -- 

%3 . Of  m 1 , . . . 


%3 . Of- 
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cerlrec,  shotnum, srcloc, sourcesize*l .25/2 . 2 , source type ,  srcz, min (dist ) , max 
(dist) ) ; 

else  %not  C-4  shot 
plottitle  =  ... 

sprintf ( ' TX02  --  Record  %g  --  %s  %s  %3.1f  m  --  %3.0f-%3.0f  m',... 
shotnum, srcloc, sourcetype, srcz, min (dist) ,max (dist) ) ; 
end 

ax=axes ('Units',  ' Normal ',  'Position',  [0.075  0.075  0.85  0.85],... 

'Visible' , 'of ' ) ; 

set (get (ax, 'Title ' ) , 'Visible ' , ' on ' ) ; 

title  (plottitle,  ' Font Size ' , 18 ,  ' FontWeight ' ,  ' Bold ' ) ; 

set (gcf , ' PaperUnits ' , ' inches ' ) ; 

set (gcf, ' Paper Posit ion ', [0.5,0.25,10,7]); 

set (gcf, ' PaperOrientation ' , ' landscape ' ) ; 

return 

tx02tbcorr — Finds  true  shot  instant 

function  [ itimeTB, itimeO , itimel , timeO , xdist , xtime]  = 
tx02tbcorr (shotnum, . . . 
cerlrec, xx, deltat, delayms, srcloc) 

%  d:\DataWinNT\FY02LoneStarTX\Plot\tx02tbcorr.m 
%  Finds  TB  (time  break  or  shot  instant)  from  data 
%  This  is  done  by  finding  the  slope  of  the  travel  time 
%  curve  for  the  first  and  last  pressure  sensor  in  the 
%  array  to  find  the  acoustic  velocity.  This  velocity 
%  is  then  used  to  find  the  shot  time. 

%  The  data  are  shifted  so  that  time=0  is  the 
%  actual  shot  time. 

%  Writes  info  to  Timebreaklogtx02 . m 
%  Does  not  call  any  other  subroutines 

%  itimeTB  is  the  TB  (shot  instant)  index 
%  itimeO  is  the  index  to  start  the  save  (usually=TB) 

%  if  the  shot  was  at  TCI,  1  second  is  added  to  itimeO 
%  itimel  is  the  index  2  seconds  after  the  TB 
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%  timeO  is  the  start  time  for  saving  the  time  series 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 


f idl3=f open ( ' Tx02TimeBreakLog . m 1 ,  1  a  1 ) ; 

srate  =  1/  deltat;%  number  of  samples  per  second 

%  Pressure  channels  to  process 
[a  b]  =  max (xx ( : , [4  8  12  16  20] ) ) ; 
b2  =  b  +  delayms*srate/1000 ;  b2  =  b2/srate; 


xp  —  [  30  60  90  120  150] ; 
if  ( srcloc (3) == 1 1 1 ) ; 

xp  =  xp  +  445  -  30;  %TC1 
elseif  (srcloc (3) =-f 2 1 ) 
xp  =  xp  +  91  -  30;  %TC2 
elseif (srcloc (3 ) == 1 3 1 ) 
xp  =  xp  +  0;  %TC3 
elseif (srcloc (3) ==*  4 1 ) 
xp  =  [281  251  221  191  161];  %TC4 
else 

xp  —  [  30  60  90  120  150] ; 
end 

vel  =  (xp  ( 5 ) -xp ( 1 ) )  /  (b2(5)  -  b2  (1) )  ; 
tlactual  =  xp(l)/vel; 
tOcorr  =  tlactual  -  b2(l); 

b3  =  b2  +  tOcorr;  %Corrected  travel  times 

%  Code  for  pistol  shots  only  -  uses  different  sensors 
if (shotnum  <=  19) ;  %Blank  pistol  shots 
if (shotnum  <=4)  %pistol  at  30  m 
iiil  =  1;  iii2  =  3;  %start  chan  4,  end  chan  12 
elseif (shotnum>4  &  shotnum<9)  %pistol  at  90  m 
iiil  =  3;  iii2  =  5;  %start  chan  12,  end  chan  20 
elseif (shotnum>9  &  shotnum<14)  %pistol  at  60  m 
iiil  =  2;  iii2  =  4;  %start  chan  8,  end  chan  16 
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elseif (shotnum>14  )  %pistol  at  30  m 

iiil  =  1;  iii2  =  3;  %start  chan  4,  end  chan  12 

end 

vel  =  (xp (iii2) -xp (iiil) )  /  (b2(iii2)  -  b2 (iiil) ) ; 

tlactual  =  2/vel;  %Pistol  shot  was  2  m  from  nearest  sensor 

tOcorr  =  tlactual  -  b2 (iiil) ; 

b3  =  b2  +  tOcorr;  %Corrected  travel  times 

fprintf (fidl3, . . . 

'Shot  %g  %s  Blank  pistol  Velocity  =  %6.1f  m/ s  TB  corr  =  %6.4f  \n',... 
shotnum, srcloc, vel ,  tOcorr )  ; 
else  %C-4  shot 

fprintf (fidl3, . . .  %C-4  shot  print  statement 

'Shot  %g  % s  CERL  %s  Velocity  =  %6.1f  m/s  TB  corr  =  %6.4f  \n',... 
shotnum, srcloc, cerlrec, vel, tOcorr) ; 
end  %code  for  pistol  shots 

itimeTB  =  -f loor ( (tOcorr  +  delayms/1000 )  *  srate); 

itimeO  =  -(tOcorr  +  delayms/1000)  *  srate; 
itimeO  =  floor ( itimeO ) ; 

timeO  =  0.0;%  Start  ascii  data  at  0.0  sec  (shot  instant) 

%  For  TCI,  start  ascii  data  at  1.0  sec 
if  (srcloc  ==  'TCI'); 
itimeO  =  itimeO  +  1.0  *srate; 
timeO  =  1.0; 
end 

if (shotnum==39  |  shotnum==4 4 ) ;  %Bad  records,  missed  shot 
tOcorr  =  0;  %Use  original  zero  time 
itimeTB  =  -f loor ( (tOcorr  +  delayms/1000)  *  srate); 
itimeO  =  -(tOcorr  +  delayms/1000)  *  srate; 
itimeO  =  floor ( itimeO ) ; 

timeO  =  0.0;%  Start  ascii  data  at  0.0  sec  (shot  instant) 
end 

itimel  =  itimeO  +  2*srate;%  Write  2  seconds  of  ascii  data 

xdist  =  xp; 
xtime  =  b3; 


f close ( f id!3 ) ; 
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return 

tx02writenchan.m  -  Writes  multichannel  ASCII  file  tx02nnn.asc 

function  tx02writenchan  (shotnum,  cerlrec,  xx,  npts,  deltat,  nchan,  delayms, ... 
descalingf act ,  .  .  . 

stackcount ,  acqdate ,  acqtime ,  engf act , xgain, itimeO , itimel ,  timeO ,  .  .  . 
dist, sensor, rtype, sensorx, sensory, sensorz, sensor2, . . . 
srcloc, stype, sourcetype, sourcesize, srcx, srcy, srcz, source2) ; 

%  d:\DataWinNT\FY02LoneStarTX\Plot\tx02writenchan.m 
%  Writes  nchan  (24)  channel  ascii  file  for  Texas  2002  data 
%  File  name  is  tx02xxx.asc,  where  xxx  is  a  new  file  number 
%  xxx  =  CERL  Shot  number  if  it  was  a  C-4  shot 
%  Calls  no  subroutines 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

srate  =  1/  deltat;  %  number  of  samples  per  second 
recnum  =  shotnum;  npts  =  length (xx) ; 
stype  =  1;  signal=l;  srccomment=source2 ; 
if (shotnum<20) ;  stype  =  5;  signal  =  0;  end 
comment  =  'Lone  Star  AAP,  Texarkana,  TX  Blast  Noise  Through  Forest 
w/CERL' ; 

%  Construct  ascii  file  name  and  open 
cerlrec2  =  str2num (cerlrec) +2000 ; 
if  (cerlrec~=0) ; 

fname2  =  sprintf ( ' tx0%g . asc ' , cerlrec2 ) ;  %tx02001  -  tx02037 
else 

if (recnumClO) 

fname2  =  sprintf (' tx02spistol0%s . asc ', num2str (recnum) ) ;  %tx02pistol0 1  - 
09 

else 

fname2  =  sprintf (' tx02pistol%s . asc ', num2str (recnum) ) ;  %tx02pistoll0  - 
19 
end 
end 

fid5=fopen (fname2, ' w' ) ; 
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%  (itimeO : itimel )  timeO  =  start  time  in  sec 

%  File  header  -  Same  for  all  records  in  this  expt : 

projectnum  =  200202;  project  =  'TX02'; 

itbtype  =  1;  tbtype= ' Manual '; %  manual; 

itbcorr  =  1;  %TB  corrected?  l=Yes 

npts2  =  16000;  %Npts  to  write  to  ascii  file 

%  File  header  -  for  all  channels  in  this  record  or  shot 
tbcorr  =  itimeO;  %Find  TB  correction  from  data  -  from  tbcorr2 
starttimes  =  timeO;  %Start  time  in  sec,  for  TCI  -  from  tbcorr2 
%stype  =  1;  signal=l ; srccomment= ' In  Open  Field'; 

%  Write  each  channel 

[aa  bb]  =  max (xx) ;  %set  plot  windows  because  manual  TB 
[cc  dd]  =  min (xx) ;  %set  plot  windows  because  manual  TB 
%for  iii  =l:nplot; 

%  Loop  over  sensors  ********************************************* 
for  iii  =  1 :nchan  % 

%  File  header  -  for  this  channel  only 
ichan  =  iii;  irtype  =  2;  reccomment  =  sensor2 ( ichan, : ) ; 
tmax  =  delayms/1000  +  (bb(iii)  ) /srate;  %times  in  seconds 
tmin  =  delayms/1000  +  (dd(iii)  ) /srate;  %times  in  seconds 
xxmax  =  aa(iii);  xxmin  =  cc(iii); 

%  Write  header 
fprintf (fid5, . . . 

'Project  no  Filename  Recno  Cerlrecno  Date  Time  \n'); 

fprintf ( fid5 ,  ' %g  %s  %s  %g  %s  %s  %s  \n',... 

pro j  ectnum, project, f name 2 , recnum, cerlrec, acqdate, acqtime) ; 

fprintf (fid5,  ' TBtype  Corrected?  Delayms  OrigNpts  TBtype  \n'); 
fprintf (fid5, . . . 

' %g  %g  %g  %g  %s  \n',... 

itbtype, itbcorr, delayms , npts , tbtype) ; 

fprintf (fid5,  'Nchan  Npts2  Srate  TBcorr  Starttime (s)  \n'); 

fprintf ( fid5 ,  '  %g  %g  %g  %g  %7.4f  \n',... 

nchan, npts 2 , srate, tbcorr, starttimes ) ; 
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fprintf (fid5,  ' Stype  Signal?  Src  (x,y,z)  Size  \n'); 
fprintf (fid5, . . . 

'  %g  %g  % 7 . 2 f  %7.2f  %7.2f  %g  %s  %s  %s  \nf,... 
stype, signal, srcx, srcy, srcz, sourcesize, srcloc, . . . 
sourcetype, srccomment) ; 

fprintf ( fid5 ,' Ichan  Rtype  Rec  (x,y,z)  \n'); 

fprintf (fid5, . . . 

'  %g  %g  % 7 . 2 f  %7.2f  %7.2f  %s  %s  \n',... 

ichan, rtype (ichan) , sensorx (ichan) , sensory (ichan) , . . . 

sensorz (ichan) , sensor (ichan, : ) , reccomment) ; 

fprintf (fid5, . . . 

'Dist(m)  Engf act , Xgain, Stack  Pmax,  T  Pmin,  T(s)  \n'); 
fprintf (fid5, . . . 

1 %7 . If  %14 . 6e  %14 . 6e  %g  %11.3e  %7.4f  %11.3e  %7.4f  \n',... 
dist  (ichan) , engfact (ichan) , xgain (ichan)  ,  .  .  . 
stackcount (ichan) , xxmax, tmax, xxmin, tmin) ; 

fprintf ( fid5 ,' %s  \n ', comment ) ; 

%  Write  data 

%  x  is  raw  data  (mV),  xx  is  data  in  physical  units, 
fprintf (fid5, . . . 

' %15 . 7e  %15 . 7e  %15.7e  %15.7e  %15.7e  \n',  ... 

xx ( itimeO : itimel , ichan) ) ; 

nwrite  =  ( itimeO+itimel+1 ) /5 ;  nwriteint  =  floor (nwrite) ; 
if (nwrite-nwriteint~=0 ) ;  fprintf ( f id5 , '  \n '  ) ; end 

end  %  loop  over  sensors  ************************ 

f close ( f id5 ) ; 

return  % - 

readnz.m — Reads  NZ  seismograph  binary  data  file 

function  [ scan, samplesPerScan, samplinglnterval , nbOf Traces ,  ... 

delaytime, descalingf act , stackcount , acqdate, acqtime]  =  readnz  (filename) 

%  Reads  a  binary  seismograph  file  in  SEG-2  format,  and  returns  the 
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%  data  in  "scan"  variable  above. 

%  This  code  written  by  Don  Albert,  CRREL,  and  based  on  code  "SEG2L0AD" 
written 

%  by  Piece  PY  to  read  a  radar  file  in  modified  SEG-2  format.  Some  of 
Piece  1 s 

%  original  code  is  still  here  in  this  file. 

%  Calls  subroutine  words. m 

%  Writes  info  about  the  file  to  file  readnzlog.m 

%  To  read  one  NZ  file  directly: 

%  fname  =  '10077.dat'; 

%  [x,  npts ,  del  tat,  nchan,  npts2 ,  delayms ,  descalingf  act ,  acqdate, ... 

%  acqtime]  =  readnz (fname) ; 

%  READNZ  -  For  Geometries  NZ  seismograph 

%This  version  modified  from  MOUT  data  for  MN02  data 

%  READNZ  Modified  by  D  Albert  to  read  standard  SEG-2  format 
%  SEG2L0AD  Read  a  SEG-2  (standard  SEG-2  format  of  the  Society  of 
%  Exploration  Geophysicist)  file  from  disk. 

%  [scan, samplesPerScan, samplinglnterval, shaf tlnterval, 

%  timerFrequency]  =seg21oad  ('filename')  reads  the  file  'filename' 

%  and  returns  the  image  scan  [m,n]  containing  n  A-Scan  of  m  samples. 

%  If  no  extension  is  given  for  the  filename,  the  extension 
%  ' . sg2 '  is  assumed. 

%  samplesPerScan  contains  the  number  of  samples  per  A-scan 
%  shaf tlnterval  contains  the  distance  between  shaft  encoder 
%  triggers  in  meter 

%  samplinglnterval  contains  the  time  between  2  samples  in  pico-seconds 
%  timerFrequency  contains  the  frequency  of  A-scan  sampling  in  Hz 

o 

o 

%  Piece  PY  24/07/1996 

%  LAMI  -  DeTeC  Demining  Technology  Center 
%  Swiss  Federal  Institute  of  Technology  (EPFL)  - 
%  Lausanne,  Switzerland 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

%  check  argument  and  filename 
if  (nargin==0) 
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error  ( ' readseg2  requires  at  least  a  filename  as  an  argument  !'); 
end; 

if  (isstr  ( filename) ~=1 ) 

error  ('Argument  is  not  a  filename  !'); 

end; 

if  (isempty  (findstr  (f ilename, ' . ' ) ) ==1 ) 

filename  =[ filename, '. sg2 '] ; 

end; 

f id=f open  (filename, ' rb ' , ' ieee-le ' ) ; 
if  (fid  ==-l) 

error  (['Error  opening  ', filename,'  for  input  !']); 
end; 

%  check  for  SEG-2  file  type 

%  first  2  bytes  equal  '3a55h'  (14933)  for  PC/Windows 

f ileType=f read  ( fid, 1 ,' short ') ; 
if  (fileType  ~=  14933) 
fclose  (fid) ; 

error  ('Not  a  SEG-2  file  !'); 
end; 

%  Open  a  log  file  to  keep  track  of  files  that  were  read 
fid2  =  f open  (  ' readnzlog . m ' ,  ' a ' )  ; 

%  fprintf ( f id2  ,  '  \n'); 

samplesPerScan  =0; 
shaf tlnterval  =0; 
samplinglnterval=0 ; 
timerFrequency  =  0; 

%  read  File  Descriptor  Block 

revNumber  =  fread  ( fid, 1 ,' short ') ; 

sizeOf TracePointer  =  fread  ( f id,  1 ,  ' ushort ' ) ; 

nbOfTraces  =  fread  ( fid, 1 ,' ushort ') ; 

sizeOfST  =  fread  ( f id,  1 ,  ' uchar ' ) ; 

firstST  =  fread  ( fid, 1 ,' char ') ; 

secondST  =  fread  ( fid, 1 ,' char ') ; 

sizeOfLT  =  fread  ( fid, 1 ,' uchar ') ; 

firstLT  =  fread  (fid,  1 ,  ' char ' ) ; 

secondLT  =  fread  ( fid, 1 ,' char ') ; 
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reserved  =  fread  ( f id, 18 , ' uchar 1 ) ; 
tracePointers  =  fread  ( fid, nbOf Traces ,  ’ ulong ') ; 

%  read  free  strings 

f seek  ( f id, 32+sizeOf TracePo inter , '  bof ' ) ; 
offset  =  fread  (fid, 1, 'ushort 1 ) ; 

%  File  descriptor  block 
while  (offset  >  0) 

freeString  =  setstr  (fread  (fid, of f set-2 ,  1  char  1 ))  1 ;  % 
%Decode  file  descriptor  block 

if  (findstr  ( freeString, ' ACQUIS ITION_DATE 1 )  >  0) 

acqdate  =  (freeString  . . . 

(length  (' ACQUIS ITION_DATE  '  ) : length  (freeString))); 
end 

if  (findstr  ( freeString, ' ACQUIS ITION_TIME 1 )  >  0) 

acqtime  =  (freeString  . . . 

(length  (' ACQUIS ITION_TIME  '  ) : length  (freeString))); 
end 

offset  =  fread  ( fid, 1 ,  ' ushort '  ) ; 
end; 

%  read  traces 

%find  number  of  samples  per  trace 
%  First  trace  descriptor  block 
fseek  ( fid, tracePointers  (l),'bof'); 
traceld  =  fread  ( fid, 1 ,' ushort ') ; 
sizeOfBlock  =  fread  ( fid, 1 ,' ushort ') ; 
sizeOfData  =  fread  ( fid, 1 ,' ulong ') ; 
nbOfSamples  =  fread  ( fid, 1 ,' ulong ') ; 
samplesPerScan  =  nbOfSamples; 

dataCode  =  fread  ( fid, 1 ,' uchar ') ; 
reserved  =  fread  (fid, 19,  ' uchar ')  ; 
offset  =  fread  ( fid, 1 ,' ushort ’) ; 
while  (offset  >  0) 

freeString  =  setstr  (fread  ( fid, of f set-2 ,' char '))' ; 
if  (findstr  ( freeString, ' SAMPLE_INTERVAL ' )  >  0) 

samplinglnterval  =  str2num  (freeString  . . . 
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(length  ( ' SAMPLE_INTERVAL  ' ) : length  (freeString) ) ) ; 
end 

if  (findstr  ( freeString, ' DELAY ' )  >  0) 

delaytime  =  1000*str2num  (freeString  . . . 

(length  ('DELAY  '): length  (freeString))); 
end 

if  (findstr  ( freeString, ' DESCALING_FACTOR ' )  >  0) 
descalingf actor  =  str2num  (freeString  . . . 

(length  ( ' DESCALING_FACTOR  '): length  (freeString))); 
end 

offset  =  fread  (fid, 1, 'ushort 1 ) ; 
end; 

%Write  file  descriptor  to  log  file 
fprintf ( f id2 , ' %s  %s  %s  %g  %g  %g  %g  \n',... 
filename, acqdate, acqtime, nbOf Samples , . . . 
samplinglnterval , delaytime, descalingf actor ) ; 

scan  =  zeros  (nbOf Samples , nbOf Traces ) ; 
for  i=l : nbOf Traces , 

fseek  ( f id, tracePointers  (i),'bof!); 
traceld  =  fread  (fid, 1, 'ushort ') ; 
sizeOfBlock  =  fread  (f id, 1 ,' ushort ') ; 
sizeOfData  =  fread  (f id, 1 ,' ulong ') ; 
nbOfSamples  =  fread  (f id, 1 ,' ulong ')  ; 
dataCode  =  fread  (f id, 1 ,  1 uchar ') ; 
reserved  =  fread  (fid, 19, 1 uchar ') ; 

%  Trace  descriptor  blocks 
%Reads  all  at  once 

freeString  =  setstr  (fread  ( fid, sizeOfBlock-32 ,' char '))' ; 

w  =  words (freeString) ; 

%MOUT  data 

%descalingf act ( i )  =  str2num (w (4, 1 : 14)  )  ; 

%stackcount (i)  =  str2num (w (18, 1 : 2 ) ) ; 

% [aasizew  bbsizew]  =  size (w) ; 

%if  (aasizew==30 ) ; 

%  stackcount ( i )  =  str2num (w (27, 1 : 2) ) ; 

%  e  1  s  e 

%  stackcount ( i )  =  str2num (w (2 6, 1 : 2 ) ) ; 

o 
o 


end 
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%stackcount ( i )  =  str2num (w (2 6, 1 : 2 ) ) ; 
descalingf act ( i )  =  str2num (w ( 7 , 1 : 14 ) ) ; 
stackcount (i)  =  X ; 

%  read  the  data  here 

scan  (1 : nbOf Samples , i)  =  tread  ( fid, nbOf Samples ,  ' float32 '  ) ; 
end; 

fclose  (fid) ;  %close  data  file 

return 

words. m — Separates  string  into  individual  words 

function  all_words  =  words ( input_string) 

%  filename:  words. m 
%  Separates  words  in  a  long  string 
%  Individual  words  are  then  in  all_words ( 9 , : ) 

%  From  the  book.  Using  Matlab,  p.  11-12 

remainder  =  input_string; 
all_words  =  ' ' ; 

while  (any (remainder ) ) 

[chopped, remainder ]  =  strtok ( remainder ) ; 
all_words  =  strvcat (all_words , chopped) ; 
end 


return 
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APPENDIX  C:  MATLAB  PROGRAMS  TO  READ  AND  PLOT 
MULTICHANNEL  ASCII  DATA  FILES. 

This  section  provides  a  listing  of  the  MATLAB  programs  used  to  read  and 
plot  the  24  channel  ASCII  data  files  that  were  constructed  from  the  binary  files 
using  the  programs  in  the  previous  section.  The  files  that  are  listed  in  this  section 
are: 

9.  doplotasc.m  -  Main  program  to  read  multichannel  ASCII  file  and  make 
a  three-panel  plot  of  the  data.  To  use,  set  the  desired  ASCII  file  numbers 
in  the  variable  “recs”  and  run  the  program.  You  may  have  to  change  the 
lines  that  construct  the  filenames  to  use  the  correct  directory  for  your 
computer,  and  you  may  want  to  change  the  print  options  at  the  very  end 
of  this  program  to  match  your  printer  or  delete  the  plot  from  the  screen 
after  plotting.  Calls  the  programs  listed  below. 

10.  readascii.m  -  Function  to  read  in  the  data  from  the  multichannel  ASCII 
files.  Calls  function  words.m 

1 1 .  plotasc.m  -  Function  to  make  a  three-panel  plot  of  the  data  in  landscape 
orientation.  The  channels  plotted  in  each  of  the  panels  are  given  by  the 
contents  of  iplotl,  iplot2,  and  iplot3. 

12.  words.m  -  Function  to  parse  a  string  into  individual  words.  This  func¬ 
tion  is  listed  in  the  previous  section. 

doplotasc.m — Reads  and  plots  multichannel  ascii  files 

%  d:\DataWinNT\FY02LoneStarTX\Plot\doplotasc.m 
%  Reads  multichannel  ascii  file,  makes  3  panel  plots  of  data 
%  Calls  functions 

%  readascii.m  -  to  read  the  multichannel  ascii  data 
%  words.m  -  to  parse  ascii  header  lines 
%  plotasc.m  -  make  a  3  panel  plot  of  the  data 

%  Ascii  record  numbers; 

%  All  C-4  (1,5,6,10,18  bad)  and  all  blank  pistol 

%  recs  =  [11  12  13  14  103  105  110  115  ];  %Example  files 
recs  =  [  2:4  7:9  11:17  19:37  43:48  101:119]; 

%  Loop  to  read  and  plot  the  files 
for  i  =  1 : length ( recs ) ; 
recnum  =  recs (i)  ; 
shotnum  =  recs (i)  ; 
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recnum  =  recnum+2000 ; 

fname  =  sprintf (' txO%g . asc ' ,  recnum) ;  %tx02001  -  tx02037 
%  Read  in  the  data  from  the  ascii  multichannel  file 
[pro j ectnum, pro j ect, fname2 , recnum, cerlrec, acqdate, acqtime, . . . 
i tbtype, itbcorr ,  delayms , npts ,  tbtype ,  . . . 
nchan, npts2 ,  srate, deltat, tbcorr, starttime, . . . 

stype, signal, srcx, srcy, srcz, sourcesize, srcloc, sourcetype, srccomment, . . . 
rtype, sensorx, sensory, sensorz, sensor, sensor2, . . . 
dist, engfact, xgain, stackcount, xxmax, tmax, xxmin, tmin, xx] =  ... 
readascii (recnum, fname) ; 

%  3  panel  plot  of  ascii  data 

%  Assign  the  channels  to  each  panel  of  the  plot 
iplotl  =  [ 9  10  11  12  5  6  7  8  1  2  3  4 ] ;  %Channels  for  1st  (left)  panel 

iplot2  =  [22  21  23  24  18  17  19  20  14  13  15  16  ] ;  %  2nd  (center)  panel 

iplot3  =  [24  20  16  11  12  7  8  3  4];  %Channels  for  3rd  (right)  panel 

%  Make  the  plot 

plotasc (shotnum, fname, xx, iplotl , iplot2 , iplot3 , deltat, starttime, . . . 
dist, sensor, sensorx, sensory, sensorz,  .  .  . 
srcloc, sourcesize, sourcetype, srcx, srcy, srcz) ; 

%  Print  and  delete 
%  print  -dljetplus; 

print  -dljetplus;  pause (10);  h=gcf;  delete (h) ; 

%  print  to  PostScript  file  for  conversion  to  PDF: 

%  print  -dps  -rl50  -append  tx02datafigs 
%  h=gcf; delete (h) ; 

end 

return  % - 

readascii. m — Reads  in  data  from  multichannel  ASCII  file 

function  [pro j ectnum, pro j ect , fname2 , recnum, cerlrec, acqdate, acqtime, . . . 
i tbtype, itbcorr , delayms , npts , tbtype, . . . 
nchan, npts2 , srate, deltat , tbcorr, starttime, . . . 
stype, signal, srcx, srcy, srcz, sourcesize, srcloc, sourcetype, srccomment, . . . 
rtype, sensorx, sensory, sensorz, sensor, sensor2, . . . 
dist, engfact, xgain, stackcount, xxmax, tmax, xxmin, tmin, xx] =  ... 
readascii (recnum, fname) ; 
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%  d:\DataWinNT\FY02LoneStarTX\Plot\readascii.m 
%  Reads  in  ascii  seismograph  data  from  fid6 
%  Calls  words. m  -  parse  string  into  individual  words 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

%  Read  first  header  to  get  number  of  channels 
f id6=f open ( fname, ' r 1 ) ; 

linel  =  f gets ( f id6 ) ; %Pro j ect  no  Filename  Recno  Cerlrecno  Date  Time 
line2  =  fgets (fid6) ; 

line3  =  f gets ( f id6 ) ; %TBtype  Corrected?  Delayms  OrigNpts  TBtype 
line4  =  fgets (fid6); 

line5  =  f gets ( f id6 ) ; %Nchan  Npts2  Srate  TBcorr  Starttime (s) 
line6  =  fgets (fid6) ; 

line7  =  f gets ( f id6 ) ; %Stype  Signal?  Src  (x,y,z)  Size 
line8  =  fgets (fid6); 

line9  =  f gets ( f id6 ) ; %Ichan  Rtype  Rec  (x,y,z) 
linelO  =  fgets (fid6); 

linell  =  f gets ( f id6 ) ; %Dist (m)  Engf act , Xgain,  Stack  Pmax,  T  Pmin,  T(s) 
linel2  =  fgets (fid6); 
linel3  =  f gets ( f id6 ) ; %Comment 
a6  =  words (line6)  ; 

nchan  =  str2num (a6 (1 , : ) ) ;  npts2  =  str2num (a6 (2 ,  : ) )  ; 

fid6=f open ( fname, 1 r '  ) ;  %rewind  file 

for  iii  =  1 :nchan  %  Read  data  ******************* 

%  Read  header  -  Same  for  all  records  in  this  expt : 

linel  =  fgets (fid6) ;  line2  =  fgets (fid6) ;  line3  =  fgets (fid6); 

line4  =  fgets (fid6) ;  line5  =  fgets (fid6) ;  line6  =  fgets (fid6); 

line7  =  fgets (fid6) ;  line8  =  fgets (fid6) ;  line9  =  fgets (fid6); 

linelO  =  fgets (fid6) ;  linell  =  fgets (fid6) ;  linel2  =  fgets (fid6) ; 
linel3  =  fgets (fid6) ;  %Comment 

%  decode  header 

a2  =  words ( line2 ) ;  a4  =  words ( line4 ) ;  a6  =  words ( line6 ) ; 
a8  =  words ( line8 ) ;  alO  =  words ( linelO ) ;  al2  =  words ( linel2 ) ; 
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%pro j  ectnum, project, f name 2 , recnum, cerlrec, acqdate, acqtime 
projectnum  =  a2  (1,  : ) ;  project  =  a2  (2,  : ) ;  fname2  =  a2  (3,  : ) ; 
recnum  =  str2num (a2 ( 4 , : ) ) ;  cerlrec  =  str2num (a2 (5, : ) ) ; 
acqdate  =  a2(6,:);  acqtime  =  a2(7,:); 
shotnum  =  recnum; 

%i tbtype, itbcorr , delayms , npts , tbtype 

itbtype  =  str2num (a4 ( 1 , : ) ) ;  itbcorr  =  str2num (a4 (2 , : ) ) ; 

delayms  =  str2num (a4 (3, : ) ) ; 

npts  =  str2num (a4 ( 4 ,  :  )  )  ;  tbtype  =  a4(5,:); 

%nchan, npts 2 , srate, tbcorr , start times 

nchan  =  str2num (a6 ( 1 , : ) ) ;  npts2  =  str2num (a6 (2 , : ) ) ; 
srate  =  str2num (a6 ( 3 , : ) ) ; 

tbcorr  =  str2num (a6 ( 4 , : ) ) ;  starttime  =  str2num (a6 (5 , : ) ) ; 
deltat  =  1/srate; 

%stype, signal, srcx, srcy, srcz, sourcesize, srcloc, 

%  sourcetype, srccomment 

stype  =  str2num (a8 ( 1 , : ) ) ;  signal  =  str2num (a8 (2 , : ) ) ; 
srcx  =  str2num (a8 (3, : ) ) ; 

srcy  =  str2num (a8 ( 4 , : ) ) ;  srcz  =  str2num (a8 (5 , : ) ) ; 

sourcesize  =  str2num (a8 ( 6, : ) ) ; 

srcloc  =  a8(7,:);  srcloc  =  words (srcloc) ; 

sourcetype  =  a8(8,:);  sourcetype  =  words (sourcetype) ; 

srccomment  =  a8(9,:); 

%ichan, rtype (ichan) , sensorx (ichan) , sensory (ichan) , 

%  sensorz (ichan) , sensor (ichan, : ) , reccomment 
ichan  =  str2num (alO (1, : ) )  ; 

rtype (ichan)  =  str2num (alO (2 , : ) ) ; 
sensorx ( ichan)  =  str2num (alO (3 , : ) ) ; 
sensory ( ichan)  =  str2num (alO ( 4 ,  : ) )  ; 
sensorz (ichan)  =  str2num(al0(5,:)); 

sensor  ( ichan, 1 : length (alO ) )  =  alO (6, 1 : length (alO) ) ; 
reccomment  =  alO (7,  :  )  ; 

sensor2 ( ichan, 1 : length (alO ) )  =  alO (7, 1 : length (alO) ) ; 

%dist (ichan) , engfact (ichan) , xgain (ichan) , stackcount (ichan) , 

%  xxmax , tmax , xxmin , tmin 

dist  (ichan)  =  str2num (a!2 ( 1 ,  : ) )  ; 
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engf act ( ichan)  =  str2num (al2 (2 ,  : ) ) ; 
xgain(ichan)  =  str2num (al2 (3,  : )  )  ; 
stackcount ( ichan)  =  str2num (al2 ( 4 , : ) ) ; 

xxmax (ichan)  =  str2num (al2 (5, : ) ) ;  tmax (ichan)  =  str2num (al2 ( 6 ,  : ) ) ; 
xxmin (ichan)  =  str2num (al2 ( 7 , : ) ) ;  tmin (ichan)  =  str2num (al2 ( 8 ,  : ) ) ; 

%  Put  data  in  xx 
data  =  [  ]  ; 

data  =  [data;  fscanf (fid6,  f  %g 1 ) ] ; 
xx (:, ichan)  =  data; 

end  %  loop  iii  =  1 :nchan  %  Read  data  ******************* 
f close ( f id6 ) ; 

return  % - 

plotasc.m — Three-panel  plot 

function 

plotasc (shotnum, fname, xx, iplotl , iplot2 , iplot3 , deltat, starttime, . . . 
dist, sensor, sensorx, sensory,  sensorz,  .  .  . 
srcloc, sourcesize, sourcetype, srcx,  srcy,  srcz)  ; 

%  d:\DataWinNT\FY02LoneStarTX\Plot\plotasc.m 

o 

o 

%  Plots  nchan  channels  of  ascii  data 
%  Does  not  call  any  subroutines 

o 

o 

%  Don  Albert  dalbert@crrel.usace.army.mil 
%  USA  ERDC-CRREL,  72  Lyme  Road,  Hanover,  NH  03755 

nplotl  =  length ( iplotl ); nplot2  =  length ( iplot2 ); nplot3  =  length ( iplot3 ) ; 
xxmax  =  max (abs (xx) ) ;  %abs  max  for  plotting 

srate  =  1/  deltat ; %number  of  samples  per  second 

%  Find  time  window  to  plot 

if  (srcloc  (3)  ==  '4')  %TC4  opposite  side  of  sensors 
tstartpanell  =  starttime  +  0.5;  tendpanell  =  tstartpanell  +  0.6; 
tstartpanel2  =  starttime  +  0.4;  tendpanel2  =  tstartpanel2  +  0.6; 
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tstartpanel3  =  tstartpanel2 ;  tendpanel3  =  tendpanel2; 
elseif  (srcloc  (3 )  ==  '1')  %TC1 

tstartpanell  =  starttime  +  0.2;  tendpanell  =  tstartpanell  +  0.6; 

tstartpanel2  =  starttime  +  0.4;  tendpanel2  =  tstartpanel2  +  0.6; 

tstartpanel3  =  tstartpanell;  tendpanel3  =  tendpanell; 
elseif  (srcloc  (3 )  ==  '2')  %TC2 

tstartpanell  =  starttime  +  0.2;  tendpanell  =  tstartpanell  +  0.6; 

tstartpanel2  =  starttime  +  0.4;  tendpanel2  =  tstartpanel2  +  0.6; 

tstartpanel3  =  tstartpanell;  tendpanel3  =  tendpanell; 
elseif  (srcloc  (3 )  ==  '3')  %TC3 

tstartpanell  =  starttime  +  0.0;  tendpanell  =  tstartpanell  +  0.6; 

tstartpanel2  =  starttime  +  0.2;  tendpanel2  =  tstartpanel2  +  0.6; 

tstartpanel3  =  tstartpanell;  tendpanel3  =  tendpanell; 
else  %Pistol 

tstartpanell  =  starttime  ;  tendpanell  =  tstartpanell  +  0.6; 
tstartpanel2  =  starttime  ;  tendpanel2  =  tstartpanel2  +  0.6; 
tstartpanel3  =  tstartpanell;  tendpanel3  =  tendpanell; 
end 

%  time  series 

tend  =  length (xx (:, 1 )) ;  %all  data 
t  =  starttime  +  ( 1 : tend) /srate ; 

%  Start  of  data  plot  %*************************************** 

nametext  =  [fNZ  Data  -  3  panel'];  % - FIGURE  -  3  Panels 

fl  =  figure (' Name  1 , nametext ) ; 

set (gcf , ' Units ' , ' inches ' , ' PaperUnits ' ,  ' inches ' ) ; 

subplot (1,3,1) ;  %Left  (1st)  panel  of  plot  - 

iplot=iplotl ;  nplot=nplotl ; 

%iplot=l:12;  %Channels  1-12  %iplot=12 : -1 : 1 ;  %Channels  1-12 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanell-0 . 05  tendpanell  -1  2*nplot] ; 
axis  (ax) ; 

axl  =  gca;  set (axl ,  ' YTick ' ,  [ ] ) ;  set (axl, 'Box' , 'on' ) ; 
xlabel ( ' Time,  sec') 

%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

plot  (t,  (xx (1 : length (t) , iplot  (i) ) /xx max (iplot (i) ) )  +y shift (i)  ,  ' k' ) ; 
leftlabel  =  sprintf (' %g  \  round (dist (iplot (i) ))) ; 
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%  Text  labels  for  plot 
strl  =  sensor ( iplot ( i ),:) ; 

str2  =  sprintf ( ’ %7 . 2e ' ,  (xxmax (iplot (i) ) ) ) ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( ' %4 . lfE-6 ' , (IE 6* xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 
str2=sprintf ( ' %4 . 2 f E-3  '  ,  ( lE3*xxmax ( iplot ( i ) ) ) )  ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( f  %4 . lfE-3 ' ,  (IE 3* xxmax (iplot (i) ) ) )  ; 

elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( ' %4 . 2  f 1 ,  (xxmax ( iplot ( i ) ) ) ) ; 

elseif  (xxmax ( iplot ( i )) )  >  1000.0 

str2=sprintf ( 1 %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 

else 

str2=sprintf ( ' %4 . If ' , (xxmax (iplot (i) ) ) ) ; 

end  %text  loop 

str3  =  char (strl); 

str4  =  char(str3  ,  str2); 

str5  =  cellstr (str4 ) ; 

text ( tendpanell  +  0 . 02,  [y shift (i) ] Astr5,  ' Font Size  1 , 8,  .  .  . 

' HorizontalAlignment 1 ,  ' left 1 ) ; 

end  %plot  loop  i=l:nplot  %Left  (1st)  panel  of  plot  - 

subplot  ( 1 , 3 , 2 ) ;  %Center  (2nd)  panel  of  plot - 

iplot=iplot2 ;  nplot=nplot2 ; 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanel2-0 . 05  tendpanel2  -1  2^nplot] ; 
axis (ax) ; 

axl  =  gca;  set  (axl ,  ' YTick 1 ,  [ ] ) ;  set  (axl ,  ' Box ' ,  ' on ' ) ; 
xlabel ( ' Time,  sec ' ) 

%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

plot(t, (xx (1 : length (t) , iplot (i) ) /xxmax (iplot (i) ) )  +y shift (i) , ' k' ) ; 
leftlabel  =  sprintf (’ %g round (dist (iplot (i) ))) ; 

%  Text  labels  for  plot 
strl  =  sensor ( iplot ( i ),:) ; 

str2  =  sprintf ( 1 %7 . 2e 1 ,  (xxmax ( iplot ( i ) ) ) )  ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( 1 %4 . lfE-6 ' , (lE6^xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 
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str2=sprintf ( ' %4 . 2  fE-3 ' ,  ( IE 3* xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf (  '  %4 . lfE-3 ' , (lE3*xxmax (iplot (i) ) ) ) ; 

elseif  (xxmax ( iplot ( i )) )  <  10.0 

str2=sprintf ( ' %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 

elseif  (xxmax ( iplot ( i )) )  >  1000.0 

str2=sprintf ( 1 %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 

else 

str2=sprintf ( f  %4 . If  1 ,  (xxmax (iplot (i) ) ) ) ; 

end  %text  loop 

str3  =  char(strl); 

str4  =  char(str3  ,  str2); 

str5  =  cellstr ( str4 )  ; 

text ( tendpanel2+0 . 02,  [y shift (i) ] ,  str5, ' Fonts ize ',8, . . . 

1 HorizontalAlignment 1 , ' left ' ) ; 

end  %plot  loop  i=l:nplot  %Center  (2nd)  panel  of  plot  - 

subplot  (1^  3 ,  3) ;  %Right  (3rd)  panel  of  plot - 

%  Pressure  sensors  only 
iplot=iplot3 ;  nplot=nplot3 ; 
yshift  =  (l:nplot)*2  -  2; 

ax  =  [ tstartpanel3-0 . 05  tendpanel3  -1  2*nplot] ; 
axis (ax) ; 

axl  =  gca;  set (axl ,  1 YTick 1 ,  [ ] ) ;  set  (axl,  'Box' ,  'on' ) ; 
xlabel ( 1  Time,  sec  1 ) 

%  plot  data  and  label 

hold  on; 

for  i=l : nplot 

if (xxmax ( iplot ( i ) ==0 ) ) ; xxmax ( iplot ( i ) ) =1 ; end 

plot  (t,  (xx (1 : length (t) , iplot (i) ) /xxmax (iplot (i) ) )  +y shift (i)  ,  ’ kf  )  ; 
leftlabel  =  sprintf (’ %g round (dist (iplot (i) ))) ; 

%  Text  labels  for  plot 
strl  =  sensor ( iplot ( i ),:) ; 

str2  =  sprintf ( 1 %7 . 2e  f ,  (xxmax ( iplot ( i ) ) ) ) ; 
if  (xxmax ( iplot ( i )) )  <  0.0001 

str2=sprintf ( ' %4 . lfE-6 ' , (lE6*xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.001 
str2=sprintf ( 1 %4 . 2  fE-3 1 ,  ( IE 3^ xxmax ( iplot ( i ) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  0.1 

str2=sprintf ( 1 %4 . lfE-3 ' , (lE3^xxmax (iplot (i) ) ) ) ; 
elseif  (xxmax ( iplot ( i )) )  <  10.0 
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str2=sprintf ( ' %4 . 2  f ' ,  (xxmax ( iplot ( i ) ) ) ) ; 

elseif  (xxmax ( iplot ( i )) )  >  1000.0 

str2=sprintf ( ' %4 . If k 1 , (  xxmax (iplot (i))/1000  )); 

else 

str2=sprintf ( 1 %4 . If ' , (xxmax (iplot (i) ) ) ) ; 

end  %text  loop 

str3  =  char(strl); 

str4  =  char(str3  ,  str2); 

str5  =  cellstr ( str4  )  ; 

text ( tendpanel3  +  0 .02,  [y shift (i) ] , str5,  1  Font Size  1 , 8,  .  .  . 

1 HorizontalAlignment 1 , ' left ' ) ; 

%  Write  title  for  right  panel 
if  (i==nplot) 

text ( tstartpanel3+0 .05, y shift (i) +1 . 1, ' Pressure  Sensors', . . . 

' VerticalAlignment ' , ' bottom ' ) ; 
end 

end  %plot  loop  i=l:nplot  %Right  (3rd)  panel  of  plot  - 

%  Title  for  entire  plot  - 

if (sourcesize~=0 )  %C-4  shot 
plottitle  =  ... 

sprintf ( ' TX02  --  %s  Rec  %g  --  %s  %3.1f  kg  %s  %3.1f  m  agl  --  %3.0f-%3.0f 

m'  ,  .  .  . 

fname, shotnum, srcloc, sourcesize*l .25/2 .2, source type, srcz , min  (dist ) , max  (d 
ist)  )  ; 

else  %not  C-4  shot 
plottitle  =  ... 

sprintf (' TX02  --  Record  %g  --  %s  %s  %3.1f  m  --  %3.0f-%3.0f  m',... 
shotnum, srcloc, sourcetype, srcz, min (dist) ,max (dist) ) ; 
end 

ax=axes ('Units',  ' Normal ',  'Position',  [0.075  0.075  0.85  0.85],... 

'Visible' , 'of ' ) ; 

set (get (ax, 'Title ' ) , 'Visible ' , ' on ' ) ; 

title  (plottitle,  ' Font Size ' , 18 ,  ' FontWeight ' ,  ' Bold ' ) ; 

set  (gcf ,  ' PaperUnits ' ,  ' inches ' ) ; 

set (gcf, ' Paper Posit ion ', [0.5,0.25,10,7]); 

set (gcf, ' PaperOrientation ' , ' landscape ' ) ; 
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return 


